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A  GENERAL  ANALYTICAL  SOLUTION  TO  THE 
PROBLEM  OP  PHOTOGRAMMETRY 

ABSTRACT 

An  analytical  treatment  of  the  general  problem  of  photogrammetry  is 
developed  which  only  requires  that  object  pointy  center  of  projection  and 
corresponding  image  point  are  collinear  for  any  one  ray  present  in  a  specific 
photogramme trie  measuring  procedure.  The  corresponding  formulas,  expressing 
the  geometrical,  relation  existing  between  the  spatial  coordinates  of  the 
object,  the  plate  coordinates  of  the  corresponding  image  and  the  elements  of 
orientation,  are  the  bases  for  a  rigorous  least  squares  adjustment,  whereby 
both  the  image-coordinates  measured  on  a  comparator,  and  the  given  control 
data  may  be  considered  as  erroneous.  The  corresponding  system  of  normal 
equations,  which  is  shown  to  be  typical,  for  any  photogrammetric  measuring 
problem^ ■ is  used  to  form  a  system  of  reduced  normal  equations,  the  vinknowns 
of  which 'are  either  the  elements  of  orientation" or  the  spatial  coordinates  of 
the  model '  The  -appllcation  of  matrix  calculus  not'  only  "simplif  ie s  the-presen- 
tatiun  bul-reduceB  the  ""booJskeepln^  effort" ,  while  ■establishing"  the -fcorre- 
spondlng  program  for  electronic  computers.  A  special  chapter  deals  wi-th  the 
problem  of  incorporating  additional  geometric  conditions  as  they  may  exist 
between  any  one  or  all  of  the  unknowns  of  the  solution.  The  application  of 
"the  described  method  for  the  special  cases  of  "strip  and  block  triangulation” 
is  treated,  whereby  attention  is  given  to  the  problem  of  solving  the  corre¬ 
sponding  system  of  reduced  normal  equations.  For  an  unlimited  strip,  a 
rigorous  solution  is  presen"ted,  which  is  based  on  the  stepwise  elimination  of 
certain  groups  of  the  unknowns .  Finally,  an  iterative  solution  is  given  which 
makes  use  of  the  relaxation  method. 
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I.  UJTRODUCnOW 


Theoretical  and  practical  work  concerned  with  the  development  and  testing 
of  analytical  solutions  for  the  reduction  of  photograrametric  records  has,  in 
the  past,  proven  the  feaslhlllty  of  such  a  method.  Such  a  project  was  initi¬ 
ated  at  the  Ballistic  Research  Laboratories,  Aberdeen  Proving  Ground,  Maryland 
with  the  study  of  the  orientation  of  a  single  camera  followed  by  a  solution 
for  the  problem  of  trlangulatlon  by  a  combination  of  two  photograrametric 
cameras.  These  two  phases  have  been  presented  in  the  BRL  Reports  No.  fl8o  ^l]* 
and  No.  961  QsJ.  The  phase  logically  to  follow  is  the  development  of  ^ 
analytical  solution  for  the  n-camera  case.  The  problem  of  triangulating,  for 
example,  the  spacQ-tlme  coordinates  of  points  of  trajectories  recorded . simul¬ 
taneously  on  more  than  two  photograrametric  cameras  belongs  in  this  category, 
as  well  as  the  problems  of  strip  and  block  triangulation. 

In  pj  an  approach  was  outlined  which  would  have  made  it  possible  to 
extend  the  two  camera  solution  to  the  treatment  of  the  n-camera  problem.  The  • 
most  serious  objection  to  that  approach  is  the  fact  that,  in  such  a  case,  the 
matrix,  associated  with  the  vectors  of  the  residuals,  loses  more  and  more  its 
diagonal  character.  Thus,  the  basis  for  the  feasibility  of  the  corresponding 
numerical  solution  is  being  impaired.  (Compare  page  22  and  schematics 

on  page  48).  Furthermore,  the  solution  becomes  cumbersome  if  additional 
Information  must  be  introduced  expressing  certain  geometric  conditions  con¬ 
cerned  with  the  coordinates  of  the  points  of  the  object  to  be  triangulated. 
Last  but  not  least,  the  "bookkeeping  effort"  in  the  preparation  of  the 
electronic  computations  would  have  been  considerable,  due  to  the  necessity  of 
distinguishing  between  various  kinds  of  control  data  and  corresponding  con¬ 
ditional  equations,  which,  to  make  things  worse,  are  of  somewhat  different 
character,  depending  on  the  number  of  cameras  Involved  in  any  one  specific 
trlangulatlon. 

The  following  solution  overcomes  these  objections.  As  an  additional 
feature,  this  solution  is  based  in  its  entirety  on  the  simplest  mathematical 
presentation  conceivable.  The  xise  of  matrix  calculus  for  setting  up  a  system 
of  reduced  normal  equations  is .especially  suited  for  electronic  computers. 

*  Reference  at  the  end  of  the  paper. 


The  following  presentation  of  the  solution  is  complete  with  respect  to 
derivation  of  the  necessary  formulas,  differential  coefficients  and  auxilia¬ 
ries.  However,  with  respect  to  general  remarks  and  hasic  principles,  referring 
to  photogreamietric  as  well  as  to  statistical  suhiects,  reference  will  occasion¬ 
ally  he  made  to  [ll  and  which  are  considered  to  he  avallahle  to  the  casual 
reader. 

II.  THE  GENERAL  PROBLEM  OP  ANALYTICAL  PHOTOGRAMMETRY 

The  general  problem  of  analytical  photogrammetry  may  he  defined  as  the 
simultaneous,  restitution  of  the  orientations  of  any  number  of  photogrammetric 
records  and  the  reconstruction  of  the  object  space  by  triangulating  corresp¬ 
onding  rays.  There  must  he  no  limitations  as  to  the  type  and  orientation  of 
any  one  camera,  so  long  as  the  corresponding  bundle  of  rays  corresponds  to 
the  principle  of  a  central  perspective.  No  limitations  must  be  made  with  re¬ 
spect  to  the  number,  type  and  location  of  control  data  including  absolute 
given,  partial  given  and  relative  points,  so  long  as  the  given  information 
satisfies  at  least  the  geometrical  req^uironents  for^  a  unique  solution. 

In  a  general  solution  it  must  be  possible  to  enforce  any  number  of  geo¬ 
metric  conditions  concerned  with  any  one  or  all  orientation  parameters,  as  well 
as  with  any  one  or  all  coordinates  of  the  points  of  the  object  space.  Further¬ 
more,  it  must  be  possible  to  consider  both  the  plate  measurements  and  the  given 
control  data  as  erroneous,  whereby  the  computation  must  allow  the  introduction, 
of  individual  weighting  factors.  Finally,  a  general  solution  must  derive,  from 
the  treatment  of  the  redundant  information,  such  expressions  of  precision  which 
will  give  Information  about  the  mean  error  of  an  observation  of  unit  weight  - 
the  mean  errors  of-  the  orientation  elements  as  well  as  of  the  triangulated 
coordinates. 


CONDITION  NECESSARY  TO  SOLVE  THE  PROBLEM  - 
OF  ANALYTICAL  PHOTOGEAMMETRY 

For  later  reference,  a  description  of  a  few  fundamental  photogrammetric 
operators  seems  appropriate. 
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A.  The  photogrammetrlc  bundle  and  its  orientation 


By  connecting  any  identifiable  detail  on  the  photographic  plate  with  a 
point  situated  outside  of  the  plane  of  the  photograph,  it  is  possible  to 
construct  a  bundle  of  rays.  (Fig.  1)  The  point,  thus  common  to  all  rays,  is 
the  so  called  center  of  projection,  denoted  by  0.  It  is  hardly  possible,  and 
for  our-  purposes  cert^i^y  unnecessary,  to  Interpret  this  point  physically. 

The  position  of  thisy^ relative  to  the  photographic  plate  is  fixed  if,  for 
example,  an  arbitrarily  oriented  rectangular  plate  coordinate  system  (x  and  y) 
is  introduced  in  which  the  base  of  the  perpendicular  to  the  plane  of  the  plate 
through  0,  denoted  by  P  (principal  point),  has  the  coordinates  x^  and  y^.  The 
length  of  this  perpendicular,  denoted  by  c,  is  a  scale  factor  and  must  not  be 
physically  Interpreted.  The  sole  purpose  of  this  phase  of  the  photograrametrlc 
evaluation  technique  is  to  provide  means  suited  for  an  unambiguous  construction 
of  a  bundle  of  rays.  This  is  accomplished  by  measuring  x,  y  plate  coordinates 
of  identifiable  detail  on  the  photograph  in  the  same  plate  coordinate  system, 
in  which  the  point  P  is  being  described  by  its  parameters  and  y^.  The 
parameters  c,  x^  and  y^  are  commonly  referred  to  as  the  elements  of  Interior 
orientation. 


After  such  a  bundle  is  obtained,  the  problem  is  to  orient  it  unambiguously 
in  space,  (l)  by  assigning  to  the  point  0,  a  specific  spatial  position  ex¬ 
pressed,  for  example,  by  three  linear  parameters,  X^,  and  with  respect  to 
an  arbitrarily  established  Cartesian  spatial  coordinate  system^  (2)  by  de¬ 
fining  with  respect  to  the  axes  of  this  coordinate  system,  the  direction  of  the 
vector  as  formed  by  the  extension  of  the  line  PO  into  space  by  two  rotational 
components,  (e.g. ,  a  and  m  rotations),  and  (5)  by  establishing  the  spatial 
position  of  the  plate  coordinate  system  by  a  third  rotation  (swing  angle  K  ) 
around  the  vector  described,  and  positioned  according  to  (2). 


The  parameters  X^,  Y^,  Z^,  a,  m,  K  are  commonly  referred  to  as  the 
elements  of  exterior  orientation. 


Both  groups,  the  elements  of  interior  and  the  elements  of  exterior 
orientation,  shall  from  now  on  be  considered  Jointly  if  reference  is  made 
simply  to  the  elements  of  orientation.  The  notation  0  is  used  if  ni  i 
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orientation  elements  are  "being  considered  present  in  a  specific  photogrammetrlc 
measuring  system,  while  the  notation  0  j  is  used  if  only  a  group  associated 
with  a  certain  camera  station  I  is  "being  regarded, 

B.  The  control  data  and  the  associated  residuals. 

It  is  common  practice  to  refer  to  points  known  with  their  absolute  po¬ 
sition  as  absolute  control  points,  to  points  known  with  either  one  or  two 
positioning  parameters  as  partial  control  points,  and  to  points  not  known  at 
fl.1 1  with  respect  to  their  absolute  position,  as  relative  control  points.  Be¬ 
cause  ftTi  these  points  belong  to  the  object  photographed,  it  is  conceivable  to 
refer  to  them  in  their  totality  simply  as  control  points,  and  denote  their 
spatial  positions  by  X  ,  while  the  spatial  location  of  any  one  individual 
point  J  shall  be  denoted  by  X  j  .  In  case  the  given  control  point  coordinates ■ 
are  not  flawless,  it  will  become  necessary  in  a  general  solution  to  treat  these 
parameters  like  measurements  and  allow  for  corresponding  corrections  in  the 
least  squares  solution.  The  designation  V  refers  to  the  corrections  of  all  . 
control  data  and  correspondingly  V  j  Is  being  used  if  only  a  specific  control 
point  "j"  is  considered. 

C.  The  plate  measurements  and  the  associated  residuals. 

As  described  in  (a),  any  individiml  ray  within  a  photogrammetric  bundle 
is  unambiguously  positioned  by  two  plate  measurements  x  and  y.  In  analogue  to 
the  above  notation,  we  shall  denote  all  plate  coordinates  involved  in  a  spe¬ 
cific  photogrammetric  measuring  system  by  x  and  the  plate  coordinates  of  a 
specific  image  point, . on  a  specific  plate  hy  Xjj .  Correspondingly,  we  denote 
the  total  number  of  plate  measurements  and  their  residuals,  by  and  V 
respectively,  and  the  plate  measurements  and  their  residuals  associated  with 
a  specific  point,  on  a  specific  photograph  by  jijj  and  Vjj  . 

The  object  space  (the  model)  can  be  visualized  as  the  integrated  effect 
produced  by  the  intersections  of  corresponding  rays,  each  of  which  is  po¬ 
sitioned  In  a  specific  oriented  photogrammetric  bundle  according  to  the 
corresponding  Image  location  as  stored  on  the  specific  photographic  record. 
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With  the  notation  introduced  above  we  may  consequently  write: 

X--.^(0,x)  (1) 

Formula  (l)  expresses  the  general  problem  of  photogramme try.  Indicating  a 
functional  relation  between  the  model  X  ,  the  plate  coordinates  X  and 
the  orientation  parameters  0  •  Introducing  the  plate  measurements,  their 

residuals  and  the  residuals  of  the  given  control  data,  we  obtain,  with  respect 
to  (l),  a  system  of  observational  equations  of  the  form: 

f(v,V)  =  F(0,X,£)  (2) 

The  roots  of  the  corresponding  normal  equation  system  represent  the  numerical 
solution  of  the  general  photogramme trie  problem.  The  result,  as  expressed  by 
formula  (1),  has  been  explained  as  an  Integrated  effect  of  the  contribution 
of  all  the  individual  rays.  Because  no  one  ray  distinguishes  Itself  baslcsLlly 
from  einy  other  ray,  this  Interpretation  suggests  that  for  an  individual  ray, 
a  corresponding  functional  relation  exists  which  is  obtained  by  simply 
indexing  the  corresponding  parameters,  thus,  leading,  according  to  formula 
(2),  to  the  corresponding  observational  equations ; 

t(Vij,Vj)  =  F(0|,Xj,iij)  (3) 

To  comprehend  this  result  we  recollect  that  the  bundle  of  rays  of  an  idealized 
physical  photograiiimetrlc  camera  has  its  geometrical  representation  in  the 
concept  of  the  central  perspective .  Any  photogramraetrlc  bundle  can  thus  be 
considered  as  a  population,  the  members  of  which  are  the  Individual  rays.  Any 
algebraic  expression  representing  a  single  ray  of  such  a  bundle,  may  thus  be 
envisaged  as  representing,  collectively,  the  bundle  in  its  entirety,  by  simply 
omitting  the  index  denoting  the  specific  ray. 

Because  any  one  photogramme trie  bundle  is  based  on  the  concept  of  the 
central  perspective,  and  any  one  photograrametrlc  problem  may  be  considered  as 
a  combination  of  any  number  of  such  bundles,  it  follows  that  an  equation 
representing  the  geometrical  properties  of  an  individual  ray  can  be  considered 
as  adequate  to  express,  collectively,  the  problem  of  analytical  photogramraetry . 
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Consequently,  it  must  be  possible  to  develop  an  analytical  solution  for  the 
most  general  problem  of  photogrammetry,  based  solely  on  formulas  which  express 
the  geometrical  properties  of  an  individual  ray  belonging  to  a  bundle  of  rays, 
formed  according  to  the  concept  of  the  central  perspective. 

The  corresponding  relation  is  the  condition  that  the  center  of  projection 
0,  the  image  point  r  and  the  object  point  R  are  collinear,  (Fig«  2) 

From  Figure  2  we  obtain: 


R  =  n  ,  r,  where  is  a  scale  factor  (it) 

The  projection  of  the  vectors  "r"  and  R  respectively  into  the  three  coordinate 
planes  gives  the  component  equations ; 


X  =  X  +  lau 

o 

Y  -  +  \iy 

Z  =  Z  +  LIW 

o 


(5) 


The  triplet  of  formulas  (5)  is  the  analytical  expression  for  the  condition 
that  the  points  0,  r  and  R  lie  on  a  straight  line , 

By  eliminating  the  scale  factor  n  in  formulas  (5)  we  obtain: 

where  (x)  =  X  -  X 


u 


(X)  -  (Z)  - 


(Y)  '  (z)  z 

W 


(y)  =  Y  -  Y 


(Z)  =  Z  -  Z^ 


(6) 


From  Fig.  2  we  read  directly: 


A_ 


r  =  iu  +  jv  +  kw  =  ix  +  jy  +  kc 


(7) 
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whereby  the  traiisfonnatlon  matrix  of  the  two  vector  triplets  is: 


1 

J 

k 


A 

1 


J. 


A 

k 


li 

A 

J1 

A 

ki 


i3 

JJ 

kj 


(8) 


Denoting,  for  ex^ple,  the  two  rotational  angles  between  the  two  vector 
triplets  according  to  Figure  2  with  a  and  o),  we  obtain: 

u  =  X  cos  a  -  y  sin  a  sin  co  +  c'‘sln  a  cos  to 
V  =  y  cos  CD  +  c  sin  co 

w  =  -X  sin  a  -  y  cos  a  sin  co  +  c  cos  a  cos  cd 
From  Figure  5  we  obtain; 


(9) 


X  =  -(x  -  X  )  cos  K  -  (y  -  y  )  sin  k 
y  =  -(x  -  X  )  sin  K  +  (y  -  y  )  cos  k 


(10) 


where  x  and  y  are  the  plate  coordinates  of  eui  Image  point  in  eui  arbitrarily 
oriented  rectangular  reference  system  (plate  coordinate  system). 

Substituting  (lO)  into  (9)  and  using  (6)  we  have: 

(z)  [  (x-x^)a^  +  (y-yp)A2  +  c  d] 


X  = 


Y  = 


Q 


(z)-[(x-x^)Bj^  +  (y-yp)Bg  +  c  e] 

Q 


^  V 


+  Y 


(11) 


and 


with  Q  =  (x-Xp)C^  +  (y-yp)C2  +  c  F 

c'  [(X)A^  +  (Y)B^  +  (zyc^^] 


X  = 


y  = 


+  X 


<1  p 

,c  [(x)A2  +  (Y)B2  +  (Z)C2] 


(12) 


+  y. 


(X)  D  +  (Y)  E  +  (Z)  F 


with  q 


where  we  have  denoted  as  essentially  direction  cosines; 

=  -cos  a  cos  K  +  sin  a  sin  oi  sin  k 
=  -cos  0)  Bln  K 

=  sin  a  COB  K  +  COB  a  Bln  tu  sin  k 

Ag  =  -cos  a  Bln  k  -  sin  a  sin  cs  cob  k  (15) 

Bg  =  cob  to  cos  K 

Cg  =  sin  a  sin  k  -cob  a  sin  to  cos  k 
D  =  sin  a  cos  to 
E  =  sin  to 
F  =  cos.  a  cos  to 

Each  of  the  pair  of  formiilaa  (ll)  and  (12)  represent  the  algebraic 
'  expression  for  the  geometrical  condition  of  co- linearity  of  the  points  0,  r 
and  R,  solved  explicitly  for  the  coordinates  of  point  R  in  the  i,  J,  k  system 
and  of  point  r  in  the  k  system^  respectively.  The  symmetrical  arreingement 

of  the  formulas  is  a  direct  consequence ,  of  the  reversibility  of  any  central 
per  spec  ti ve .  ( c  ompare ,  [ll  page  8) 

In  [2]  the  algebraic  solution  was  continued  by  eliminating  such 
and  from  the  solution  which  were  not  given  as  control  data.  This  step  led, 

for  relative  control  points,  to  an  algebraic  expression  for  the  condition  of 
intersection,  and  correspondingly  for  partially  given  control  points,  to  the 
condition  of  intersection  at  either  one  or  two  given  control  coordinates.  As  a 
direct  consequence  of  that  approach,  different  types  of  conditional  equations 
were  obtained  for  the  different  types  of  control  points.  To  complicate  the 
situation  further,  the  number  of  any  particular  set  of  such  conditional  equations 
depends  on  the  number  of  camera  stations  Involved  in  any  specific  trlangulatlon 
problem. 
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TJie  most  serious  obJectlon',agalnBt  the  approach  chosen  In  [2]  results 
from  the  fact  that,  for  the  condition  of  intersection,  It  Is  necessary  to  have 
at  least  four,  and  under  certain  conditions  even  six  plate  coordinate  measure¬ 
ments  together  with  their  residuals  in  any  single  conditional  equation. 
Furthermore,  in  the  case  of  multi-camera  trlangulation,  the  same  residuals 
appear  in  more  than  one  conditional  equation,  making  a  rigorous  least  squares 
treatment  laborious,  or  in  more  complex  cases,  even  impractical. 

The  coordinates  of  the  relative  and  partial  control  points  not  given  at 
the  outset  of  the  computations  become  known  in  the  solution,  as  presented  in 
[2J  ,  only  Impllcltlly  as  functions  of  the  adjusted  plate  measurements  and 

the  orientation  elements.  The  introduction  of  any  additional  geometric  con¬ 
ditions  as  they  may  exist  for  any  one  or  all  of  these  coordinates  requires 
complex  mathematical  manipulations,  prohibitive  from  the  standpoint  of  com¬ 
puting  economy. 

These  difficulties  have  been  avoided  in  the  present  solution  by  including 
the.  elimination  of  the  unknown  coordinates  of  the  model  in  the  process  of  the 
numerical  treatment.  In  other  words.  Instead  of  eliminating  algebraically  the 
unknown  X^,  and  Z^*s,  beforehand.  In  the  system  of  formulas  (12),  these 
quantities  eire  carried  as  unknowns  together  with  the  orientation  elements  and 
are  solved  for  during  the  numerical  solution. 

All  of  the  beneficial  consequences  of  this  approach  will  become  evident 
in  the  least  squares  solution  as  described  in  the  next  chapter.  It  suffices 
to  mention  here  the  most  important  features; 

(l)  Each  point,  without  regard  to  its  character  as  absolute,  partial  or 
relative  control  point  is  being  treated  alike,  thus  giving  rise  to  two  ^ 
equations,  of  the  form  of  formulas  (l2),  for  each  camera  station  at  which  the 
point  was  recorded.  Any  given  information  for  example,  in  the  form  of  spatial 
coordinates,  is  introduced  by  simply  eliminating  the  corresponding  parameter 
corrections  from  the  sequence  of  unknowns  in  the  least  squares  solution. 
Therefore,  all  points  are  treated  with  only  one  basic  set  of  eqviatlone, 
(formulas  (12)),  which  are  explicit  in  terms  of  the  measured  plate  coordinates 
X  and  y,  respectively. 
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(2)  The  fact  that  the  coordinates  of  the  points  of  the  model,  to  be 
triangulated  are  carried  as  unknowns  in  the  solution,  makes  it  possible  to 
Introduce  readily  any  additional  existing  geometric  conditions  for  any  one 
or  all  of  these  coordinates. 

IV,  THE  LEAST  SQUARES  SOLUTION 

The  formulas  (12)  express  the  plate  coordinates  x  and  y  as  functions  of 
the  orientation  elements  and  the  spatial  coordinates  of  the  corresponding 
object  point.  For  any  one  point  J  photographed  at  a  certain  camera  station  I, 
we  may  therefore  write  In  general  terms,  according  to  formula.  (5): 


£  +  V  =  F, 

’‘ij  ^ 

[  (otyb,  >c  ,X^,Y^,Z^,c,Xp)^,  (X  +  Vj^,  Y  +  V^,  Z  +  V2)j  ] 

(U) 

£  +  Y  =  F„ 

^  (q!,U),  k  ,X^,Y^,Z^,c,y^),^,  (X  +  Vjj,  Y  +  Vy,  Z  +  ^Z^j] 

From  the  Taylor  expansion  for  the  right  hand  side  of  the  equations  (l4), 
neglecting  terms  of  second  and  higher  order,  we  obtain  the  observational 
equations : 


.dF  .  dF  .  dF  .  dF  SF  SF  SF 


ij 

aF, 


dF, 


r( 


dF. 


■^iJ 


hF 

^F  , 

+  Ac  + 

75 — ^  ^ 

Sf„ 

Sf„ 

2  q.  ,r 

3r  c  \ 

2  q 

3Z“  c  ^Z 

SF2  S 

^2  \ 

J'p  / 1 

3f, 


dF, 


Sf, 


M 


Sf, 


ij 

hr' 


(15) 
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SF. 


dF. 


^IJ 


where 


-  Ajg  =x"~Ji 


-At  =  y“  .  i 


The  V  are  the  normalized  V-residuals  according  to  the  formula; 

c 

x°  and  y°  are  computed  with  formulas  (12)  and  the  approximation  values  of 

IJ  IJ  ft  H 

the  unknowns  which  are  denoted  by  .  It  is  obvious  that  in  any  one  specific 
set  of  obseivatlonal  eq.uatlons  may  appear  either  certain  coordinate  corrections 
or  the  corresponding  residual  errors  V  depending  on  the  cheiracter  of 
the  object  point  under  consideration.  With  the  introduction  of  the  V-reslduals, 
according  to  formula  {l6)f  each  ray  has  been  assigned  a  specific  group  of  such 
residuals.  Such  an  approach  is  desirable  from  the  numerical  standpoint,  be¬ 
cause  any  correlation  is  avoided  between  the  coefficients  of  the  matrix 
associated  with  the  residuals  of  the  various  rays,  intersecting  at  a  specific 
point.  From  the  theoretical  standpoint  such  a  solution  seems  to  be  advantageous 
because  the  individual  bundles  of  rays  will  conform  to  the  pattern  of  the  con¬ 
trol  data  without  undue  restraint.  In  paragraph  (p)  of  this  chapter  the  possi¬ 
bility  will  be  discussed  of  arranging  the  least  squares  adjustment  in  such  a 
way  that  for  each  given  control  coordinate  only  one  specific  residual  V  is 
obtained,  independent  of  the  number  of  rays  intersecting  at  the  point  imder 
consideration. 


The  relation  between 
unknowns  is  given  by: 

o 

a  =  a 


o 

03  =  03 


o 

K  =  K 


the  approximation  values 


+  A  0! 
+  A  03 

+  A  »c 

+  A  X 

0 

+  A  Y 

0 

+  A  Z 


o 

c  =  c 


and  the  final  values  of  the 


+ 

+ 

+ 

+ 

4 

4- 


A  c 
Ax 


A  X 


A  Y. 


A  Z, 


(17) 
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A,  A  Direct  Solution 


Using  Biatrlx  notation,  the  system  of  observational  equations  for  an  m-ray 
solution,  according  to  formulas  (5)  or  (15),  may  be  vritten  as: 


The  A'l  are  the  coefficient  matrices  of  the  corresponding  residual  vectors 
Vj  .  In  case  the  absolutely  given  control  data  are  considered  flawles^ 

A,  =  Ags  A3  *  •  •  •  =  I 

the  unit  matrix.  The  Bj  are  the  coefficient  matrices  of  the  vectors  of  the 
corresponding  parameter  corrections  A  j  .  The  i  j  are  the  vectors  of  the 
absolute  terms  of  the  observational  equations.  We  may  rewrite  the  system  of 
observational  equations  with  obvious  notation  as  A  V  *  B  A-Ji  (19) 

with  the  weights  P  , 

Assuming  the  observations  to  be  independent  and  normally  distributed  the  most 
probable  values  of  the  unknowns  are  obtained  by  minimizing  PV  ,  where 
P  denotes  the  weight  matrix 
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An  individual  P|  matrix  designates  the  weights  of  the  plate  measurements  and 
of  the  control  data,  as  they  pertain  to  a  specific  ray.  The  introduction  of 
weighting  factors  for  the  plate  measurements  may  hecome  necessary  to  express 
varying  degrees  of  precision  associated  with  the  original  observations.  Such 
variations  may  arise  from  the  method  of  measuring,  the  varying  image  quality 
caused  by  loss  of  definition  towards  the  edges  of  the  photograph,  or  from  a 
decrease  of  accuracy  with  which  the  distortion  correction  is  known  for  image 
points  at  increasing  radial  distances  from  the  center  of  the  plate.  Various 
degrees  of  dependability  of  the  given  control  data  can  be  considered,  as  well, 
by  the  introduction  of  corresponding  weighting  factors.  Infinitely  large 
weight  with  respect  to  the  given  control  data  will  eliminate  the  V  from  the 
solution,  thixB  distributing  the  unavoidable  discrepancies  among  the  plate  co¬ 
ordinate  measurements  only;  however,  infinitely  large  weight  assigned  to  the 
plate  coordinate  measurements  will  make  the  given  control  data  absorb  all  the 
discrepancies  present  in  connection  with  a  certain  ray  under  consideration.  [  3] 
While;;. the  boundary  conditions  that  were  Just  mentioned  are  only  of  secondary 
importance,  the  choice  of  suitable  weighting  coefficients  inakes  it  possible 
to  prevent  an  undue  deformation  of  the  model  obtained  photogrammetrically  by 
possible  strain  in  the  configuration  of  the  given  control  data.  In  case  the 
original  observations  (image  point  or  control  coordinates)  are  not  independent 
from  each  other,  the  corresponding  correlations  can  be  allowed  for  by  intro¬ 
ducing  the  corresponding  correlation  coefficients  into  the  matrix.  In 
case  B-1 1  the  p-Values  are  equal,  it  is  convenient  to  consider  P  as  the  ■unit 
matrix. 


In  setting  up  the  corresponding  normal  equations,  one  has  to  take  into 
account  the  fact  that  in  the  most  general  case  each  observational  equation 
contains  more  thsui  one  measurement.  Furthermore,  in  such  a  case,  certain 
measurements  and  their  residuals  appear  in  more  than  one  obser'vational 
equation.  ■  Helmert  in  [4j  ,  (pp  215-222),  has  shown  a  direct  solution  of  the 
general  problem  of  a  least  squares  adjus'tment.  (Compare  [^2  j  ,  paragraph  IV ). 
Accordingly,  we  obtain  in  our  case  a  set  of  normal  equations  as  shown  in 
formulas  ( 21 ) ; 


AP-*  K  -  EA  =  -ii 

'  -  k  ^  0 


(21) 
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TECHNlCALLfSKAKV 
branch  no.  ] 


vhere  k  denotes  a  vector  of  unknown  Lagrange  multipliers  or  correlates.  A 
direct  solution  would  "be  obtained  by  the  inversion  of  the  system  (21 ) .  That 
la,  hovever^  hardly  practical  because  In  addition  to  the  unJoiown  parameter 
corrections,  A»s,  ve  have  for  an  m-ray  problem  2m  additional  unknown  k- values, 
a  fact  which  increases  the  number  of  unknown  parameters  beyond  the  means  of 
practical  computations. 


Next,  a  method  which  is  commonly  used  in  partitioning,  a  system  of  linear 
eq,uatlonB  shall  be  presented  which  will  be  useful  for  our  problem.  Assuming  a 
system  of  linear  equations; 


Ax  =  JL 


(22) 


or  conseq.uently 

x=  A"'  1 


(23) 


Such  a  system  can  be  partitioned  in  any  arbitrary  manner  leading  to  the  fo31ov- 
ing  sub-matrices  and  sub- vectors; 


X 

< 

+  A|2  X2  =  ft. 

(24) 

A2|X| 

+  A2  Xj  =  52 

According  to  (23) 

i  .. 

we  may  write  with  the  notation  of  (24) 

Xt 

'< 

II 

^I  =  0i“A|2  X2 

(25) 

and  therefore 

X| 

=  A7'  ft, 

-A",'  A|2  X2 

(26) 


Introducing  (26)  into  the  lower  portion  of  the  system  (24)  we  obtain; 

^21  A  j’  di  “A21  A  Aj2  Xg  +  Ag  X2  =  ft2  (27) 

(  A2“A2i  A|'  A|2  )  X2  "  ftz'Azi  A|* 


21 


which  may  be  written  in  a  form  analogous  to  (25)  as;  ' 

^2  “  ^2  ^2  with  Ag  =  Ag  “  Agi  A|  A|2 

^2  “  -^2  "  ^21  -^l 

The  computation  of  Xj  is  then  carried  out  with  formula  (26). 

Because  the  sequence  of  the  steps  in  the  process  of  partitioning  is  In 
no  way  restricted,  it  is  possible  to  write  as  in  (28)  and  (26) 

X,  =  A*'‘  X,  with  A*  =  A, -AjgAgAg,  (29) 

X|  =  ij  “  Ajg  Ag  ig 


The  method  Just  described  obviously  eliminates  one  of  the  two  groups  of 
uisknowns  as  chosen  by  the'^process  of  partitioning  and  solves  for  the  other 
group.  If  the  method  la  used  to  partition  a  system  of  normal  equationa  at 
any  point  along  its  diagonal,  it  follows  from  the  symmetry  of  auch  a  system 
that  A|2  “  Agj  .  Further,  it  can  be  shown  that  in  such  a  case  the 

matrix  Ag  in  (28)  and  correspondingly  the  matrix  A*  in  (29)  are  again  ' 
symmetrically  arranged  square  matrices. 

Because  the  matrix  AP’^A’’’  in  formula  (2l)  is  non-rsingular  in 

our  problem,  we  may  apply  the  method  of  partitioning  as  Just  described  for  the 
purpose  of  eliminating  the  k- values  from  the  original  normal  equation  system. 
The  reduced  normal  equation  system  is,  according  to  formula  (27), 


[B^(AP"'A^r'  b]a  =  B'^{AP"'a^)"'1 
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The  feaslhllity  of  this  method  of  partitioning  depends  on  the  effort  necessary 
to  invert  the  matrix.  AP-'A^.  Because  this  matrix  is  even  for  the  most 
general  case  of  our  problem^  a  sequence  of  fully  separated  symmetrically 
arranged  (2  x  2)  square  sub-matrices,  it  is  possible  to  accumulate  the  normal 
equation  system  (51)  stepwise,  as  explained  in  [2]  page  22,  Thus  we  obtain; 

.1  [  BT(  AP-'  A“)“*B]jA=l  [B'lAP-’A'^r*  jl]j  (52) 


whereby  m,  the  number  of  AP“*  at  submatrices,  equals  the  number  of  rays 
present  In  the  specific  problem.  As  already  mentioned  at  the  beginning  of  this 
paragraph,  in  case  only  the  residuals  of  the  plate  measurements  are  present  the 
Aj  matrices  are  unit  matrices  and  therefore  the  (AP“'  A"^)"*  term  in  (31) 
reduces  to  P  .  In  such  a  case  the  system  (31)  resembles  a  system  of  normal 
equations  associated  with  ohservatlonal  equations  for  Independent  indirect 
measurements.  The  final  normal  equation  system  in  such  a  case  can  he  accumu¬ 
lated  stepwise  according  to  formula  (32)  by  considering  in  each  computational 
step  a  single  observational  equation. 


After  the  vector  of  the  A  corrections  of  the  unknown  parameters  is 
computed,  we  obtain  with  the  first  group  of  equations  in  formula  (21),  the 
k -values. 


/  A  i-»— I  aT  I  /  n  A  A  \ 
K  =  V  wr  ■  \  D4i-  li  / 


and  the  residuals  V  and  V  by: 

''x 

Vy 

Vx 
Vy 


V  = 


LVzJ 


=  P”'  A^  k 


The  V-valuea  are  then  computed  with  formulas  (16), 


(33) 


(34) 


From  a  study  of  formula  (51)  It  becomes  evident  that  the  feasibility  of 
the  proposed  solution  depends  on  the  possibility  of  Inverting  a  more  or  less 
large  matrix  of  normal  equations.  In  special  applications  of  photogrammetry, 
such  as  laboratory  measurements,  ballistic  measurements,  terrestrial  appli¬ 
cations  and,  more  recently,  cadastral  surveying,  using  airborne  photography, 
the  single  model  still  plays  a  dominant  role.  The  total  number  of  unknowns 
in  such  cases  may  very  well  be  within  the  limits  which  today  can  be  reduced 
by  Inverting  the  corresponding  normal  equation  matrix.  However,  depending  on 
the  number  of  relative  control  points  considered  essential  in  order  to  satisfy 
the  requirement  for  redundant  information  and  on  the  type  of  the  electronic 
computer  available,  there  will  arise  in  many  photogramme trie  applications  the 
problem  of  treating  numbers  of  unknown  parameters  exceeding  the  computational 
capacity  for  reduction  by  direct  inversion.  For  such  cases,  it  is  desirable, 
however,  to  maintain  the  advantages  of  the  systematlsra  and  simplicity  described 
earlier  in  forming  the  observational  equations  and  the  corresponding  normal 
equations.  Consequently,  the  solution  to  our  problem  must  concern  itself  with 
methods  for  determining  the  roots  of  a  normal  equation  system  of  the  type  shown 
in  formula  (5l). 

B.  A  Solution  by  Partitioning 

In  order  to  further  reduce  the  number  of  unknowns  in  the  normal  equation, 
system  obtained  in  formula  (51),  we  split  the  B  matrix  and  the  A  vector  in 
such  a  way  that  one  group  is  associated  with  the  model  and  the  other  group  with 
the  camera  orientations.  We  denote  the  corresponding  submatrlces  and  subvectors 
by  BxiBo  and  Ax  and  Aq  ,  respectively.  With  these  notations  we  can 
present  the  system  of  normal  equations  (51)  as  follows: 


[bJ  (AP''a^)''  Bx 

Iax 

+ 

bI  (AP''a^)  '  Bo 

Ao  =  bI  (ap^'a’^I'i 

[B5(AP"'A^r'  Bx] 

Ax 

+ 

Bo(AP''A’'r'  Bo. 

Ao=  Bj(AP''A’'r'l 

2k 


The  normal  equation  system  as  shown  In  (55)  Is  typical  In  its  arrangement  for 
any  photogrammetrlc  triangulation  problem.  The  number  of  points  recorded  at 
any  one  station  and  the  number  of  stations  involved  in  a  specific  measuring 
program  will  obviously  influence  the  overall  size  of  this  matrix  system  but 
will  not  change  its  basic  character.  A  study  of  the  corresponding  matrix 
(see  Fig.  4)  shows  that  we  have  along  the  diagonal,  a  sequence  of  fully 
separated  symmetrically  arranged  square  submatrices.  The  fact  that  the  two 
types  of  submatrices  which  appear  veiry  in  size  is  less  significant  than  the 
fact  that  in  each  submatrix  of  the  first  group,  bJ  ( AP  *  ,  the 

spatial  coordinates  of  only  one  specific  object  point  (in  a  general  case  up  to 
5)  are  present,  while  in  each  submatrix  of  the  second  group  Bo(AP  ^A"^)  ^  Bq^ 
only  the  orientation  elements  of  one  specific  camera  station  (in  a  general  case 
up  to  9),  appear.  The  BxtAP~*A^)~’  Bq  and  the  Bq  (AP~*  A^  )“*  Bx 

submatrices  express  the  fact  that  a  specific  point  was  photographed  from  certain 
camera  stations. 


Viewing  our  system  of  normal  equations  with  respect  to  the  method  of 
partitioning  as  described  in  this  report  with  formulas  (22)  through  (50),  a 
suitable  point  for  partitioning  is  obviously  that  point  on  the  diagonal,  which 
separates  the  parameters  associated  with  the  model,  from  the  parameters  con¬ 
nected  with  the  camera  orientations,  as  Indicated  by  the  dotted  lines  in 
formula  (55)«  According  to  formula  (27),  we  may  write: 


{bJ  Bo-[Bo  (AP'V^  )■'  Bx][bI(AP  Vr'  Bx]  '[^(AP‘'aV  Eb]}Ao  = 

{bJ  (AP^'a^)"'  -[b5  (AP‘'aV  Bx][Bx(AP"'a4'Bx]''  [bJ (ap''a^ r']}l 


or,  with  reference  to  (32): 


I  (CoBo)iAo  =  i  (CoA); 


(57) 
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vherc 


Co|=b5  {(AP-'ATr'-[(AP-»  ATr'Bx][Bj(AP"‘A^r'Bx]  '  [Bl(AP“'ATr']}| 


and  r  =  nvimber  of  points  of  the  model. 


T*^  rkvi^T  rfT^QTnTViQ4-  'v»-i  r*  -I  TYr»/>V>T /:iVM  T.fo 

.i«li  CUi.jr  ■♦■■«  t.i  I  <-1.1111111-^  1.  1^  I  t-  4,  s.; 


■monr  + /*\  o^^Tn■^nQ^^ 

UW  --1  -■-  t  1..  ■  1  LI.1.  t-jv 


coordinate  corrections  of  the  model  or  the  corrections  of  the  orientation 
parameters  in  the  system  of  equations  (51)*  Accordingly,  ve  may  write  as  in 
formula  (3Y) 


n  n 

2  I  Cx  Bx  )jAx’“.^  ^  (58) 

j=l  1=1 


irfiere: 

CK-B;{(Ap-'An-'-[(AP-'AT)-'Bo][B5(Ap-'ATr'Bo]''[B5(AP-'Ah-']}| 


and  n  =  number  of  photogrammetric  cameras. 


In  photogrammetric  measuring  problems,  generally  spealdLng,  the  number  of 
unknown  orientation  parameters  will  be  less  than  the  number  of  unknown  coordi¬ 
nates  of  tbe  model]  therefore  the  elimination  of  Ax  as  suggested  with 
formiilas  (37)  will  generally  lead  to  the  most  economical  solution.  Block 
triangulation  vith  a  high  degree  of  sidelap  may  be  mentioned  as  an  exception. 


With  the  vector  A  q  known,  the  vector  Ax^^sy  he  computed  for  each 


^-1- 


Thus  for  example  Ax 


A  =  r  rj  t&p-i  at  D.. 


-I 


rTjaD-IaT\-I  /  ft.D  A  \ 


f3Q'l 
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with  A 


■[J:] 


the  k  and  the  corresponding  v  and  V  values  may  he  computed 


with  formulas  (55) ^  (5^)  and  (l6). 


However,  it  appeaj^s  to  be  more  advisable  to  compute  the  coordinates  of 
model  points  with  such  computational  steps,  as  outlined  in  Chapter  IV,  E.  of 
this  report.  Such  an  approach  appears  to  be  especially  suited  because  the 
corresponding  computational  means  must  be  included  in  the  solution  anyway,  in 
order  to  provide  the  first  approximation  values  for  the  coordinates  of  the 
model,  which  are  necessary  for  starting  the  initial  iteration  cycle. 


The  setting  up  of  the  observational  equations  (15)  or  (19)  requires 
computing  the  coefficients  of  the  matrix  of  the  unknown  parameter  corrections 
and  of  the  residuals.  All  coefficients  are  obtained  by  partial  differenti¬ 
ation  of  the  formulas  (12).  The  linearization  procedure  is  accomplished  by 
applying  the  Taylor  series  and  neglecting  second  and  hi^er  order  terms. 
Therefore,  an  Iterative  procedure  must  be  provided  in  the  computation,  whereby, 
the  results  of  each  cycle  are  Introduced  as  approximate  values  to  the  follow¬ 
ing-cycle.  The  iteration  is  repeated  until  the  solution  has  converged  to  a 
pre-established  accuracy  level. 


We  now  introduce  the  following  computational  auxiliaries: 


c 


axj  c 


"''la 

^  3 


=2 


(D  -  ©  E  - 
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The  coefficients  of  the  observational  equations  (12)  are  now: 
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The  q,  k^,  Dj,  E,  F  terms  are  computed  with  formulas  (12 )j, 

(15)  and  the  approximation  values  of  the  unknown  parameters. 


It  should  be  mentioned  that,  for  special  cases  where  the  orientations  of 
the  various  camera  axes  are  such  that  the  approximation  values  for  the  ro¬ 
tational  parameters  can  be  assumed  in  all  iteration  cycles  to  be  either  0, 

^  or  multiples  there-of,  the  differential  quotients  as  given  in  (4l)  reduce 
to  well  known  simple  expressions.  However,  this  fact  is  hardly  worth  con¬ 
sidering  for  a  general  solution  because  the  savings  in  computing  time  are 
Immaterial  if  high  speed  electronic  computers  are  used. 

Small  and  medium  sized  electronic  computers  will  handle  with  this  method, 
photogtammetric  measuring  systems  of  two  and  three  camera  stations  with  an 
unlimited  number  of  recorded  points.  Large  computers  are  adequate  to  solve 
the  corresponding  problem  for  five  and  six  camera  stations.  In  other  words, 
the  suggested  method  seems  feasible  to  provide  a  practical  analytical  solution 
for  photogranmjRtrlc  measuring  problems  encountered  in  terrestrial,  laboratory 
cadastral  and  ballistic  applications.  The  problems  of  strip  and  block  tri- 
angulation  which  need  to  be  studied  further  are  considered  in  a  later  chapter. 

0.*^  On  the  Use  of  Celestial  Control  Points: 

The  use  of  celestial  targets  (mostly  fixed  stars)  is  a  traditional  means 
in  geodesy  and  photogrammetry  for  establishlrg  absolute  orientation  of  certain 
bundles  of  rays  with  respect  to  earth  fixed  coordinate  systems.  Iterrestrial 
Photogrammetry  has  used  this  method,  especially  for  non-topographlcal  appli¬ 
cations.  This  technique  has  lately  become  of  interest  again'  in  connection 
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with  tracking  of  guided  missiles  and  satellites.  In  addition,  star  photography- 
must  be  considered  an  excellent  means  of  calibrating  precision  photogramme trie 
cameras,  because  the  me-thod  not  only  allows  the  determination  of  angular 
parameters,  but  the  determination  of  -the  elements  of  Interior  orientation. 
Including  distortion  coefficients;  as  well.  (See  chapter  VI  of  -this  report). 

Whatever  -the  purpose,  the  me-thod  relies  simply  on  the  fact  -that  each 
celestial  -target  point,  known  by  its  spherical  coordinates,  right  ascension 
(a)  and  declination  (B),  provides  an  absolute  control  point  situated  at 
Infinity,  Thus,  a  direction  In  space  Is  fixed  which  is  not  sensltl-ve  to 
changes  in  -the  position  of  the  cen-ter  of  projection.  Provided  that  the  date 
and  time  of  exposure,  together  with  -the  geodetic  ellipsoidal  coordinates  of 
-the  center  of  projection  (0  and  X)  are  known,  such  a  direction  can  be  expressed 
wi-th  well-known  formulas  by  a  set  of  azimuth  and  elevation  angles,  or  as  more 
often  used,  by  so  called  standard  coordinates  |  and  tj,  which  are  referred  to  a 

■teuigent  plane  on  -the  unit  sphere,  described  around  the  cen-ter  of  projection. 

< 

With  proper  orientation,  the  £  and  tj  coordinates  In  the  notation  of  this  re¬ 
port,  are  equal  to  (X),  (y)  -wl-th  (Z)  =  +1.  (Compare  formulas  (12)) 

The  corresponding  observational  equations  are  formed  exactly  like  those 

for  absolu-te  control  points,  with  the  exception  that  each  of  -the  partial 

differential  coefficients  associated  with  "the  AX,AY,AZ  parameter 

o  o  o 

corrections  blames  zero,  because  of  the  celestial  target  points  being  at 
Infinity, 

It  ~ls ’practical  to  incorporate  the  Influence -of  aB-tronomical  refraction, 
(see  chap-ter  VII  of  this  report)  In  -the  computation  of  the  s-tandard  coordl- 
na-tes,  so  that  for  -the  corrected  rays  no  further  attention  need  be  paid  to  the 
refraction  problem  during  -the  trlangulatlon  computations . 

Obviously,  -the  sun  edso  Is  such  a  celestial  target  and  consequently,  the 
aforementioned  facts  concerning  the  data  reduction  apply  to  It  as  well.  The 
use  of  -the  sun  In  aerial  photogrammetry  Is  Important  for  strip  and  block 
trlangulatlon,  and  appears  "to  be  mandatory  for  satelli-te  photogrammetry,  as 
-well  as  for  special,  attitude  determinations  of  airborne  photogramme-tric 
cameras.  The  significance  of  celestial  -targets  is  not  res-trlc-ted  to  the  fact 
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that  a  direction  is  available  which  is  not  affected  hy  accumulated  errors. 

At  least  eq^ually  important  Is  the  fact  that  because  of  the  Infinity  position 
of  celestial  targets,  the  solution  is  able  to  discriminate  between  rotational 
and  translational  parameters.  If  only  terrestrial  control  points  are  used 
these  parameters  are  closely  correlated.  (Compare  )• 

Generally  speaking,  we  can  say  that  each  recorded  sun  Image  gives  rise  to 
two  observational  equations  in  terms  of  the  three  rotational  parameters  (a,  oi, 
k)  of  the  corresponding  station.  Obviously,  it  is  necessary  that  the 
measured  sun  image  plate  coordinates,  recorded  with  a  special  sun  camera,  be 
normalized  in  such  a  way  that  they  are  compatible  with  the  metric  character¬ 
istics  of  the  aerial  camera  taking  the  corresponding  ground  photography. 
Mathematically  speaking,  it  is  necessary  to  normalize  the  corresponding 
comparator  measurements,  by  transforming  the  x,y,c  system  of  the  sun  camera, 
with  three  rotations  and  three  translations,  so  that  it  conforms  with  the  x,y, 
c  system  of  the  aerial  camera.  The  corresponding  rotational  and  translational 
parameters  are  obtained  from  preceding  camera  calibrations,  for  example  by 
taking  star  photography. 

D.  Taking  into  Consideration  Additional  Geometric  Conditions 

In  order  to  complete  a  generally  acceptable  analytical  solution  it  is  now 
necessary  to  study  the  problem "of  incorporating  certain  geometric  conditions 
such  as  may  exist  for  one  or  all  of  the  unknown  parameters,  including  both  the 
coordinates  of  the  model  and  the  orientation  elements. 

In  the  above  outlined  solution,-  a  system  of  normal  equations  exists,  at 
least  temporarily,  in  terms  of  both  the  elements  of  orientation  and  the 
coordinates  of  the  model.  The  incorporation  of  any  additional  conditions 
existing  for  any  one  or  all  of  these  unknown  parameters,  can  be  performed  by 
a  computational  technique  which  was  presented  by  Helmert  in  fk]  dealing  in 
Chapter  IV,  paragraph  2k  on  page  196ff  with  the  problem  of  indirect  measure¬ 
ments,  the  unknown  parameters  of  which  have  to  satisfy  certain  conditional 
equations . 
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Helmertr  shows  In  his  direct  solution,  which  seems  suited  for  our  problem, 
that  it  is  only  necessary  to  add  to  the  normal  equation  system  the  corre¬ 
sponding  conditional  equations  and  to  restore  the  quadratic  form  of  the  equation 
system  by  Introducing  for  each  added  conditional  equation  an  unlmown  auxiliary. 
Assuming  the  observational  eqmtlons 

Vi  =  a^^x  +  b^y  +  c^z  - 

=  a^x  +  b^y  +  c^z  -  ^2 

V,  =  a-x  +  b,y  +  c,z  - 

5  5  5  5  5 

^  =  a^x  +  b^^y  +  c^z  - 
and  the  conditional  equations 


0  =  Pq  +  +  p^y 

0  =  q  +  q^x  +  q,z 
o  1  5 

We  obtain  from  (42)  the  corresponding  normal  equations: 


(45) 


[aa]  X  + 

+ 

>> 

1 — 1 

[ac]  z  - 

[af]  =  0 

[ab]  X  + 

[bb]  y  + 

[be]  z  - 

[bi]  =  0 

(44) 

[acj  X  + 

[cb]  y  + 

[cc]  z  - 

[cf]  =  0 

This  system  is  being  enlarged  by  adding 

the  following  form: 

the  conditional  equations 

(45),  into 

[a^  X  +  [a^ 

y  +  [ac] 

z  +  Pll^j_ 

+  ^1^2  ~ 

=  0 

[a^  X  +  (b"^ 

y  +  [bcl 

z  +  Pghj^ 

,  -  [bl 

=  0 

[ac]  X  +  [cl3 

y  +  [c(3 

z  + 

+  ^5^2  ■ 

=  0  (45) 

P^x  +  ^2^ 

o 

0 

-  +  P 

•^o 

=  0 

q^^x 

+  q^z 

9 

•  +  q 

o 

=  0 

The  presence  of  zeros  in  some  of  the  diagonal  terms  must  be  taken  Into 
consideration  when  the  system  (45)  is  being  inverted. 
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All  the  possible  conditional  eqviations  which  may  exist  between  the 
unknown  parameters  of  our*  problem  can  be  divided  into  two  types.  In  the  first 
group,  all  conditions  shall  be  considered  which  introduce  additional  metric 
information,  as  for  example,  the  height  of  a  point  above  a  certain  reference 
surface,  or  the  partial  position  of  a  point  expressed  by  latitude  0  and  , 
longitude  \  or  the  sla^t  distance  between  two  relative  control  points  of  the 
model  or  the  differences  in  flying  height  between  two  successive  aerial 
photographs,  etc.  This  type  of  conditional  equation  which  shall  be  designated 
as  a  metric  condition  is  especially  Important  for  the  wide  field  of  photo- 
grammetric  applications  using  ground  control  points. 

The  use  of  photogrammetric  measuring  systems  for  the  deteSrmlnatlon  of 
trajectories  of  missiles  and  the  anticipated  use  of  these  methods  to  measure 
precisely  the  trajectories  of  satellites,  e.g.  for  obtaining  geodetic  basic 
information,  requires  the  Introduction  of  a  second  type  of  conditional  equation, 
where  the  additional  information  is  not  metric  in  an  absolute  sense  but  re¬ 
stricted  to  the  mathematical  character  of  the  trajectory.  This  type  of 
conditional  equation  shall  be  designated  as  trend  ^  condition. 

(a)  Metric  Conditional  Equations 

The  mathematical  form  of  a  specific  metric  condition^  equation  is 
Influenced  by  the  type  of  metric  information  given,  and  by  the  reference 
coordinate  system  Introduced  for  the  solution  of  the  specific  photogrammetric 
measuring  problem.  From  the  multitude  of  possible  given  metric  information, 
the  following  selection  can  be  assumed  to  be  of  special  interest  for  terrestri¬ 
al  control  points. 

absolute  points  given  by:  Latitude  (0),  longitude  (l)  and  elevation 

(h)  with  respect  to  an  ellipsoid  of  revo¬ 
lution,  the  axes  of  which  are  denoted  by 
a  and  b 

partial  points  given  by:  either  latitude  (0)  and  longitude  (x),  or 

by  elevation  (h) 

relative  points  given  by:  the  slant  distance  between  two  relative 

control  points 


For  computational  convenience  (see  Chapter  VIl),  it  is  assumed  that  a  local 
Cartesian  system  denoted  by  XYZ  is  oriented  in  such  a  way,  that  its  XY  plane 
is  tangent  to  the  reference  ellipsoid  at  the  point  of  local  origin,  and  its 
+X  axis  is  pointing  to  the  south. 


The  corresponding  metric  conditional  equations  are  for: 

(l)  An  absolute  control  point  given  by  0,  X,  H:  X  =  (0  X  H) 


y  =  F^  (0  X  H) 
Z  =  (0  XH) 


The  system  (46)  gives  the  results  of  a  coordinate  transformation  presented  in 
[t).  No  special  metric  conditional  equations  become  necessary. 

(2)  A  partial  control  point,  given  by  0  and  X: 

(2)  I^X  +  +  III^Z  .  A,  (47) 


I^X  +  +  III^Z  . 

and  correspondingly,  if  the  approximation  values  X' 
conditional  equations 
(1) 

(2) 

where,  in  formulas  (47)  and  (48); 

^0  =  +1 
0 


°  y°  Z°  satisfy  the  above 


+  lyY  +  =  0 

=  0 


1^0  - 


(48) 


IIl^  =  cos  r01  tan  0  -  sin  m.  cos  X^ 
sin  w  tan  0  +  cos  cos  X^ 
*  0  =  X*  +  iiy* 

I  =  tan  X*  sin  [ 0 j 


with  X*  =  X  _[X] 


II 


-1 


III  ^  =  tan  X*  COB  [0] 


(49) 


A 


X 


=  -  [X]  tan  X^ 
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[x],  [y]  0,  [z]  are  the  geocentric  coordinates  of  the  local  origin  given  "by 

[0],  W.  [H]  =  0.  The  geocentric  system  is  oriented  so  that  its  +X  axis  lies 

^  ^ 

in  the  meridiem  plane  of  the  local  origin.  The  X  Y  Z  coordinates  are 
Cartesian  coordinates  in  the  local  system  for  the  point  given  hy  0,  \  with  H  =  0. 


It  should  be  noted  that  all  the  coefficients  can  be  computed  during  the 
phase  of  coordinate  transformation  and  consequently  enter  as  constants  into  the 
least  squares  solution  of  the  actual  photogrammetric  triangulation  computations, 

(5)  A  partial  control  point,  given  by  H: 

Such  a  point  is  situated  on  an  ellipsoid  of  revolution,  the  axes  of  which 
are  a  —  (a  +  H)  and  b  =  (b  +  H).  In  order  to  make  the  approximation  values 
X°  Y°  Z°,  as  they  may  be  obtained  with  (57),  compatible  with  this  condition,  it 
is  suggested  that  first  a  new  Z  value  be  computed  with: 


The  corresponding  conditional  equation  is  then 

IIjjAY  +  =  0 

whereby: 

Ig  =  +  X  '1I2  “  2, 

o 

IIjj  =  hY 

III^  =  Tlj  -  X°T1^ 

and 

T)^  =  cos^  [?^]  +  sin^[0] 

.,2  =  Sln2[j«]  +  (%)2 

b 

q,  =  [X]  cos  [0]  +  (^)^  [z]  sin  [0] 

b 


(50) 

(51) 


(52) 


(55) 


* 

=  Sin  [0]  cos  [0]  {(^)^  “  l} 
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"*6  ^  '^8  =[xlsin[0]  -[z]  cos  [0]  iiQ=[X]sin[0] -(%)^  [z]cos[0]  (53! 

'l  ‘1  ^  "b 

The  r] -values  are  constants  and  best  computed  during  the  preparatory  coordinate 

transformation  computations . 


(4)  In  the  case  that  the  slant  distance  (d)  between  two  relative  control 
points  1  and  J  is  given,  we  obtain  the  corresponding  metric  conditional  equation 

d  =  [  (X^  -  +  (Y^  -  +  (Z^  -  Zj)2] 

and  correspondingly: 

(X°-XpAX.+(Y°-YpAY^+(Z°-Zp^^-(X°-XpAX.-(Y°-YpAY.-(Z°-ZpZffij  =  A  d  (54) 

where;  „  o  ^  ri  on 


A  d  = 


The  metric  conditional  equations  (48),  (51)  and  (54)  are  linear  equations  in 
terms  of  coordinate  corrections  and  have  the  form  of  equations  given  with 
formulas  (43) ,  Consequently,  these  conditional  equations  can  be  used  by 
slnply  adding  the  specific  equations  to  the  normal  equation  system  as  given  in 
formulas  (35)^  which  in  turn  have  to  be  restored  to  a  quadratic  form  according 
to  the  system  presented  with  formulas  (45).  The  resulting  system  is  solved  by 
direct  inversion,  or  solved  by  the  formulation  of  a  system  of  reduced  normal 
equations  by  partitioning  as  described  in  the  previous  paragraph.  It  is 
obvious  that  during  the  process  of  partitioning  and  forming  the  final  normal 
equation  system,  the  minimum  size  of  any  one  subdivision  in  the  computing 
process  is  determined  by  the  specific  group  of  parameter  corrections  combined 
by  one  or  more  conditional  equations.  In  order  to  obtain  a  uniform  approach 
to  the  computing  program,  it  is  possible  to  consider  all  points  of  the  model 
as  relative  points,  introducing  the  absolutely  given  control  data  as  metric 
conditional  equations.  Considering  the  given  coordinates  as  approximate  veilues, 
the  corresponding  conditional  equations  with  respect  to  the  the  three  spatial  co 
ordinates  are  AX.=bO,  AY  =0  and  A  Z  .  =  0.  These  conditional  equations, 

J  o  J  T  1  T  I 

or  any  selection  of  them,  are  added  to  the  corresponding  B  (AP”  A' )’ 
during  the  process  of  forming  the  reduced  normal  equation  system  (37)*  By  so 
doing  the  coordinate  values  Introduced  are  enforced.  If  desired,  the  corre¬ 
sponding  V  residuals  can  be  carried  in  the  A  -matrix  as  described  before. 
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,If  on  the  other  hand,  the  conditional  equations  are  rather  simple 
expressions  as  e.g.  in  formulas  and  (5l),  it  is  more  practical 

to  use  these  equations  to  eliminate  certain  unknowns  in  the  corresponding 
observational  equations.  Thus,  the  coefficients  in  the  matrix  as  given 
with  (4l)  have  to  be  replaced  by: 


(L)^  -  -  F, 

(L)^  =  -  F^ 


If  for  points  of  type  5  the  unknowns 
A  X  and  A  Y  are  eliminated 


H 

and  (j)  =  F  -D 

'■  'x  X  III„  X 

n 

It, 

f  J)  =  F  i  -D 
y  y  11%  y. 


(k)  =  f  ytI  -  e 

X  xll%  X 

(K)  =  F  “  -  E 

'  "y  y  11%  y 


If  for  points  of  type  h  the 
unknown  AZ  is  eliminated 


In  case  the  systematic  errors  present  in  the  control  data  are  rather  large, 
it  may  be  desirable  to  adjust  them  in  such  a  way  that  only  a  specific  residual 
V  is  obtained  for  any  given  control  coordinate,  independent  of  the  number  of 
rays  intersecting  at  such  a  point.  Such  a  solution  may  be  obtained  by  adding 
to  the  normal  equation  system  (2l),  corresponding  conditional  equations.  In 
general  nomenclature,  such  a  conditional  equation  is,  with  regard  to  formulas 
(15),  (16),  (34)  and  (40), 

V=  p~^(akj^  +  bkg)  J 

Thus  it  is  possible  to  establish  for  each  of  the  independent  combinations  of 
the  intersecting  rays  by  pairs,  a  conditional  equation  between  the  two  group 
of  associated  and  k^  values,  by  equalizing  the  corresponding  expressions 
(55)  separately  for  each  given  control  coordinate.  In  this  way  all  normal 
equations  associated  with  a  specific  point  will  become  interlocked  in  the 
AP"»  Al  matrix  and  the  process  of  stepwise  accumulation  of  the  reduced 
normal  equation  system  (32)  or  (3?)  will  be  rendered  much  more  cumbersome. 
Numerically  speaking,  in  such  cases  it  appears  simpler,  and  from  the  theo¬ 
retical  standpoint  sufficiently  rigorous,  to  accomplish  the  adjustment  in 
two  separate  steps.  In  the  first  step  a  least  squares  adjustment  will  be 
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perfonned  according  to  the  presented  solutlor^  with  the  restriction  however, 
that  only  such  absolute  given  control  data  are  Introduced  as  are  necessary  for 
a  unique  solution.  The  coordinates  of  the  model  thus  obtained  are  then  trans¬ 
formed  In  a  second  computational  step  by  three  translations,  three  rotations 
and  a  scaling  factor  in  such  a  way  that  the  sum  of  the  squares  of  the  residual 
distances  between  the  model  and  all  the  given  control  coordinates  becomes  a 
minimum.  Obviously,  with  such  a  method,  the  photogramme tr leal  ly  obtained  model 
Is  being  interpolated  Into  the  configurations  of  the  given  control  coordinates. 
(Compare  [jQ  and  [6j  ) 

(b)  Trend  Conditional  Equations 

Quite  similar  to  the  treatment  of  a  metric  conditional  equation,  it 
id  possible  to  introduce  a  certain  trend  condition  known  to  exist  between  the 
points  to  be  triangulated,  and  which  Is  mathematically  expressed  hy  a  functional 
relation  existing  between  the  coordinates  of  such  points.  In  a  general  sense, 
such  conditional  equations  will  resemble  metric  conditional  equations,  except 
that  they  win  not  necessaxlly  have  any  absolutely  given  parameters.  Depending 
on  the  mathematical  character  of  a  specific  trend,  the  number  of  points  Involved 
In  such  a  trend  condition  will  vary.  Correspondingly,  the  smallest  possible 
subdivision  in  the  process  of  establishing  the  final  normal  equation  system  by 
partitioning  will  be  determined  by  the  number  of  points  combined  by  any  one  or 
several  trend  conditional  equations. 

E.  A  Solution  for  Triangulating  Points  not  Included  In  the  Least  Squares 
Treatment  for  the  Orientation  Parameters 

;  Despite  the  possibility  of  Incorporating  any  number  of  relative  control 
points  In  the  general  analytical  solution  for  a  specific  photpgrammetrlc 
trlangulatlon  problem.  It  sametlmes  may  be  desirable  to  triangulate  separately 
additional  points  of  the  model.  Consequently,  for  those  points,  an  independent 
coordinate  determination  becomes  necessary.  The  positions  of  the  corresponding 
rays  are  determined  by  the  elements  of  orientation  as  obtained  from  an  Indei- 
pendent  preceding  least  squares  adjustment  and  the  corresponding  plate  measure¬ 
ments.  It  is  obvious  that  these  rays  will  not  intersect  due  to  unavoidable 
measuring  errors.  They  must  be  to  intersect  so  that  the  sum  of  the  squares 

of  the  corrections  to  be  applied  to  the  original  piate  measurements  Is  a  minimum. 
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It  appears  advantageous  to  have  an  analytical  solution  which  is  affected  as-- 
little  as  possible  by  the  number  .of  rays  involved  in  a  specific  triangulation 
case,  and  which  in  addition  uses  formulas  already  in  use  in  the  computational 
procedures  previously  described.  Obviously,  such  a  solution  is  already  available 
with  the  observational  equations  (15)  if  only  the  coordinate  corrections  A  X  , 

J 

A  Yj  and  A  are  considered  as  unknowns .  These  observational  equations  lead 

directly  to  the  corresponding  normal  equation  system  formed  according  to 

formulas  (55 )  by  Introducing  (AP‘*A^  as  P  and  Bq  as  a  null  matrix.  The 

necessary  coefficients  of  the  observational  equations  denoted  by  J  ,  J  .  K  , 

X  y  •  X 

K  and  L  and  L  are  given  in  formulas  (4l). 

y  X  y 

This  approach  makes  it  necessary  to  compute,  as  a  first  step  approximation, 
values  for  the  coordinates  of  the  point  under  consideration.  This  may  be  done 
efficiently  with  the  use  of  the  formulas  (U),  which  may  be  writ-ten  as: 

X  +  a  Z  +  p  =  0 

X  ^  N 

(56) 

Y+aZ+S=0 

“■»r 


where ; 


(x-x  )A^  +  (y-y jAg  +  cD 


(x-x  )Bi  +  (y-y + 


p  =  -  (a  z  t  X  ) 

X  X  o  o 

p  =  -  (a  z  f  Y  ) 

'  y  y  o  o 

The  corresponding  normal  equation  system  for  an  n-ray  solution  is; 


+  [“x]  ^ 


H  -"[“yl  [Py]  -  °  (57) 

+  [oa]  +  [ap]  =  0 

It  should  be  pointed  out  that  the  roots  obtained  from  formula  (57)  must  not 
be  considered  as  the  result  of  a  rigorous  least  squares  solution  because  there 
is  no  Indication  how  nearly  this  approach  minimizes  the  sum  of  the  squares  of 
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the  corrections  of  the  original  plate  measurements.  Without  doubt  the  answer 
will  provide  an  excellent  approximation  value.  ,  As  a  matter  of  fact,  such  an 
answer  may  be  considered  as  adiequate  in  itself,,  if  there  is  evidence  of  the 
presence  of  systematic  errors,;  a  situation  which  renders  an  additional 
treatment  by  a  rigorous  least  squares  solution  superfluous.  Furthermore,  the 
method  presented  with  the  formulas  (56)  and  (57)  is  well  suited  for  computing 
in  each  Iteration  cycle  of  the  least  squares  solution  as  described  in  Chapter 
IV,  the  coordinates  of  the  points  of  the  model.  (Compare  remark  Chapter  IV  (b), 
page  no.  (28). 

P.  The  Determination  of  the  Mean  Errors  of  an  Observation  of  Unit  Weight 
of  the  Elements  of  the  Orientation  and  of  the  Triangulation  Results 

The  mean  error  of  an  observation  of  unit  weight  denoted  by  m  is  computed 

with 

„  .  (VlPV)l/2 
'  r-u 

The  term  V^Pv  ,  may  he  obtained  directly  from  the  reduction  of  the  normal, 
equations  or  by  adding  the  8q\iares  of  the  individual;  weighted  v  and  V  values . 
The  letter  r,  denotes  the  number  of  observational  equations  and  u  denotes  the 
mnnbftr  of  unknown  parameters.  Thus  the  mean  error  of  a  specific  observation  I 
before  adjustment  is; 

(59) 

The  computation  of  m  directly  from  the  original  measurements,  e.g.  using  the 
differences  of  multiple  observations,  may  lead  to  a  value  of  greater  physical 
significance.  The  discrepancies  between  the  different  values  of  m,  computed 
with  different  methods  provide  means  to  investigate  the  presence  of  systematic 
errors . 

The  mean  errors  of  the  unknown  parameters  in  a  least  squares  solution  are 
obtained  by  multiplying  m  with  the  corresponding  weighting  factors.  The 
inverse  of  the  matrix  of  the  coefficients  of  the  final  normal  equation  system, 
is  the  matrix  of  the  weighting  coefficients.  The  diagonal  elements  Eire  the 
squares  of  the  weighting  factors  associated  with  the  corresponding  unknown 
parameters . 


m,  = 


m 


If  the  normal  equation  system  as  given  with  formula  (Jl)  is  directly 
lnvei*ted,  the  computation  of  the  corresponding  mean  errors  for  both  the  orien¬ 
tation  elements  and  the  coordinates  of  the  model  does  not  cause  any  difficulty. 
Correspondingly,  the  mean  errors  of  the  unknown  parameters  associated  with  the 
reduced  normal  equation  system  in  the  method  of  partitioning  (xg  In  formula 
(28)  or  in  formula  (29),  respectively)  are  directly  obtained  using  the 
diagonal  terms  of  the  inverse  of  the  corresponding  normal  equation  matrix.  To 
prove  this,  we  write  the  matrix  of  the  coefficients  of  the  normal  equation 
system  with,  a  notation  corresponding  to  formula  (24)  as: 

(60) 


The  corresponding  Inverse  matrix  shall  be  denoted  by: 

(61) 


According  to  the  definition  of  the  process  of  Inversion,  with  respect  to  the 
notation  used  in  formulas  (60)  and  (6I)  follows: 

^12  ^11  Q|2  _  Fi  0 

A|2  A22  Q|2  ^22  ^ 
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In  accordance  with  the  definition  of  the  matrix  product,  it  follows  from  (62) : 


(1) 

A,,  Q„  + 

A|2  0^2  =  ^ 

(2) 

All  Q|2  + 

A 1  2  Q22  ■  Q 

(5) 

A12Q11  + 

A22  QT2  ■  0 

(4) 

^12 

^22  ^22 '  ^ 

With  (65)  ve  obtain: 

Q,2  =”A|]  A|2  Q22 

T  T  -I 

^12  Q|2  ■  "  ^12  ^11  ^12^22 
■H  A22Q22*  ^22  ^22 

A|2Q|2‘*‘  ^22^22  *^6^)4  =  I  -  (  A22”A[2  Al'l  A, 2)  Q22 

or 

Q22'  ^  ^22"’  ^12  ^11  ^12^  * 


(63) 


(64) 


Formula  (64)  is  identical  with  the  corresponding  expression  A  2  formula 

(28),  which  had  to  he  proven. 


In  order  to  compute  the  mean  errors  of  those  parameter  corrections,  which 
are  eliminated  during  the  process  of  forming  reduced  normal  equations,  ac¬ 
cording  to  formulas  (57)  or  (38)>  the  matrices  Q j|  or  022»  respectively,  must 
he  obtained.  It  follows  from  formulas  (65) ^  equations  (l)  and  (2),  that 

(65) 


0||  '  *'il  + 


A||'  A|2  Q 


.T  .-I 
22  ^12  ^11 
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(65) 


and  correspondingly  from  the  eqiiatjLona  (3)  and  (4) 

^22  *■  -^22  '*’^22  ^[2  ^11^12  ^Z2 


The  computation  of  the  matrices  Q||  or  Q22can  not  he  performed  hy  stepwise 
accunnilatlon  and  therefore  a  considerable  computing  effort  becomes  necessE^« 

On  the  other  hand,  the  relative  accuracies  within  a  model  are  often  of  pro¬ 
nounced  interest.  This  fact  underlines  the  advantages  of  a  solution  which  Is 
established  either  directly  on  formula  (52)  or  on  the  reduced  normal  equation 
system  (58)^  provided  that  the  corresponding  computations  can  be  handled  by 
electronic  computers.  If  It  suffices  to  obtain  only  the  mean  errors  of  the 
eliminated  parameter  corrections.  It  Is  possible  to  compute  the  squares  of  the 
corresponding  weighting  coefficients  as  the  diagonal  terms  in  (65)  or  (66)  by 
a  stepwise  accumulation  In  the  same  way  as  the  reduced  normal  equation  system 
was  stepwise  accumulated.  .  ,  “ 

In  case  the  elements  of  orientation  can  be  assumed  as  flawless,  the  corre¬ 
sponding  Qj^iDis-trix  becomes  a  null  matrix  and  the  yx'i*  matrix  emerges  as  the 
weight  matrix  of  the  triangulated  points ;  The  more  excess  observatians  in¬ 
corporated.  into  the  original  least  squares  solution  for  the  orientation 
parameters,  the  better  this  approximation  solution  will  be. 

G.  An  Example  of  the  Described  Solution. Using  Functional  Schematics 

Figure  3  shows  the  overlap  of  3  photographs  which  may  be  considered  as 
being  taken  either  by  aerial  or  ground  established  cameras.  The  different 
types  of  control  points  ere  marked  by  the  following  symbols; 

^  ;  absolute  control  point,  given  by  X,  Y,  Z 

d) :  partial  control  point,  given  by  X  and  Y 

^  j  partial  control  point,  given  by  Z 

0;  relative  point  -- 
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(66) 


The  following  tabulation  shows  the  situation  Interpreted  in  terms  of  the  re¬ 


cording  of  the  various  points  (Nos.  1-9)  by  the  various  camera 

stations  (Nos. 

I,  II 

and  III) 

Point 

lype  of  Point 

Camera  Stations 

at  which  the  point  was  recorded 

Nos. 

I 

II 

.  Ill 

1 

A 

X 

-  • 

X 

2 

0 

X 

- 

X 

3 

@ 

X 

X 

X 

4 

0 

- 

X 

X 

5 

0 

X 

-  X 

- 

6 

0 

X 

X 

X 

7 

A 

- 

X 

X 

8 

X 

X 

■  - 

9 

A 

- 

X 

- 

counting  the  cross  marks^  we  see  that  19  individual  rays  are  present.  Each 
ray  leading  to  2  equations,  we  have,  19  x  2  =  -38  observation  equations  for  the 
33  unknowns  which  (assuming  for  example,  the  three  elements  of  interior  orien¬ 
tation  of  each  camera  orientation  as  known)  are  composed  of  3  *  6  =  18  elements 
of  orientation  and  (9  x  5)  12  =15  coordinates  of  the  model.  As  an  inter¬ 

esting  by-product,  it  may  be  mentioned,  that  the  evaluation  of  the  chosen 
problem, ;  although  slightly  overdetermlned,  would  not  be  possible  with  con¬ 
ventional  restitution  equipment  and  techniques,  because  there  are  only  a 
maximum  of  4  intersections  between  corresponding  rays  of  any  one  pair  of  photo¬ 
graphs,  as  can  be  seen  from  the  above  tabulation. 

With  the  help  of  this  tabulation  we  form  the  corresponding  38  observation** 

al  equations  which  are  shown  schematically  in  Fig.  6.  In  Fig.  6  each  line 

represents  an  observational  equation  formed  according  to  formulas  (18)  with  the 

terms  given  in  formulas  (I3),  (40)  and  (4l),  The  individual  A  matrices  and 

Blmllarly  the  corresponding  P -matrices  actually  are  a  sequence  of  completely 

separated  sub-matrices.  An  Inspection  of  Fig.  6  shows  the  fact,  which  has 

« 

already  been  mentioned,  that  in  case  there  are  no  V*s,  the  A  matrices  become 
a  sequence  of  imlt  matrices. 
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/[  corrections  of  point  coordinates 

/l. 


A Q  corrections 


of  orientation  elements 


No.  of  unknown  parumefer  correction  s  ^ 


f^ig^  1 
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Absolute  terms 


The  actual  character  of  a  specific  control  point  with  respect  to  its 
absolute,  partial  or  relative  value  is  of  no  concern  in  setting  up  the  corre¬ 
sponding  observational  equations .  In  case  any  control  coordinates  are  given, 
it  Is  only  necessary  to  eliminate  the  corresponding  parameter  corrections  in 
the  Bx  -matrix.  Therefore,  no  coordinate  corrections  appear  for  the  absolute 
control  points  Nos.  1,  7  and  the  partial  control  point  5  appears  only  with 
one,  and  the  partial  control  point  8  with  only  two  coordinate  corrections. 

The  system  of  corresponding  normal  equations  formed  according  to  formula 
(31)  is  presented  in  Figure  7. 

V.  A  REFERENCE  FOR  AW  AUXILIARY 
COORDINATE  TRANSFORMATION 

The  subject  of  a  specific  coordinate  transformation  is  not  necessarily 
connected  with  the  subject  of  this  report.  However,  the  application  of  photo- 
grarametry,  especially  for  geodetic  purposes,  unavoidably  confronts  the  user 
with  the  problem  of  converting  geodetic  ellipsoidal  coordinates  (latitude 
longitude  X,  and  elevation  H,  as  referred  to  an  ellipsoid  of  revolution)  into 
a  system  of  arbitrarily  oriented  Cartesian  coordinates  and  vice  versa. 

A  solution  for  this  problem  is  given  in  [7]  under  the  title  "Some  Remarks 
on  the  Problem  of  Transforming  Geodetic  EUlpsoldal  Coordinates  into  Cartesian 
Coordinates  with  the  Help  of  the  Reduced  Latitude".  Obviously,  such  a  solution 
includes  the  establishment  of  a  geocentric  Cartesian  system,  which  is  but  a 
special  case  of  the  general  coordinate  transformation  problem. 

VI.  DETERMINATION.  OF  RADIAL  DISTORTION 

Distortion  A  is  positive  if  the  image  point  is  displaced  away  from  the 
principal  point. 


Figure  8 
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Figure  9 


From  Fig,  8: 


=  [< 


/  .2  /  .?n  1/2 


Assuming  that  the  distortion  A  can  be  expressed  by 


and 


Ax  =  ^  (x  -  X  )  =  fK  +  K  d^  4-  K  4.  TT 

a  r  p-  '  o  "1  “2~  '  “5"  ' 


U„  ..  \ 

•  •  •  <■  -  Ap/ 


(67) 

(68) 

(69) 


A5r  -  ^  (y^  -  y^)  =  (k^  +  K^d^  +  K^d^  +  Kjd^+ . )(y^  .  y^) 

Due  to  the  character  of  lens  distortion,  it  is  immaterial  if  d  is  interpreted 
as  the  radial  distance  with  or  without  distortion. 

From  Fig.  9: 


A  -  A>  ^  d 
Ac  c 


(70) 


or 


A*  =  A  -  —  d 
c 

+  ^  +T1+^  nrF  J  .«  /rr/NN 

VA  usAw.i.xj.^  Av^iuiua.cv  XII  \{^) 


A»  =  K^d  +  K^d^  +  K^d^  +  K^d'  +  .. 
wherebv 

V  . 

K'  >  K  - 

O  O  C 


(71) 


is  zero  if  the  condition  is  Introduced  that 

K  =  — 
o  c 


(72) 
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Therefore,  assuming  that  In  the  orientation  calibration  the  principal  distance 
c  Is  considered  as  unknown,  that  Is  to  say  that  the  solution  Is  not  deprived 
of  the  possibility  of  computing  a  A:-correctlon  In  each  Iteration,  we  may 
present  the  remaining  distortion  by; 


and  correspondingly  with  formulas  (69)  and  (73)  and  Pig.  8: 

\  ^  \  +  (K^d^  +  Kgd^  +  Kjd^  +  ...)(x^  -  x^) 

"  ^r  -  Yp) 


From  formulas  (12)  and  with 
and  formulas  (74)  the  observation  equations  are; 


=  Jt  +  V 


y„  =  i  +  V 

R  y  y 


(73) 


(74) 


(75) 


“  F  ""p  “  ^1 


where 


i  +y  =  ~  (1  +  K,d‘=  +  K-d^  +  K,d  +  ..,)+ 
y  y  q  1  ^  5 


c  /  2  .  2.1/2 


P 

y.,  =Fc 


(76) 


d  =  —  (m  +  n  ) 

q 

m  =  (X)A^  +  (Y)B^  +  (Z)C^ 
n  =  (X,)A2v+:(t.)B2  +  (Z)C2 

Considering  Kj^,  and  as  a  sufficient  number  of  distortion  parameters 
we  obtain  from  the  Taylor  expansion  for  the  right-hand  side  of  the  above 
equation  according  to  fonmilas  (15 )^  neglecting  terms  of  second  and  higher  order, 


"^x  ••• 

=  ^  Aa  + 

• 

V  +  b  Vv  +  . . . 

y  X  A 

=  ■5^ 

^1 


^F  dFp  SPp 


(77) 


^  ^1  ^  “‘2 
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Because  the  unknowns  K^,  Kg,  and  are  linear  in  the  observational ' 
equations  (76),  their  approximations  in  each  iteration  may  be  t^en  equal  to 
zero;  consequently  the  coefficients  of  the  remaining  unknown  orientation 
parameters  remain  the  same  as  given  by  formula  ( 41 ) . 

The  incorporation  of  the  distortion  determination  requires  merely  the 
addition  of  the  partial  differentials  with  respect  to  the  unknowns  Kj^,  Kg  and 
Kj  to  the  system  of  observational  eqxiations. 

They  are^  for  the  x- equation: 


for  K,  :  M  =  d^ 

1  X 


cm  ,2  ,  o  Ov 
—  =  d  .(£  -  X  ) 

q  '  X  p' 


=  d" 


«2  =  \  = 


K,  ;  0  =  d 

and  for.  the  y-eqmtlon: 


cm  -H-  /  .u  o\ 

—  =  d  .  (i  -  X  ) 
q  X  p 


.6 


cm  ,b  /  .o  0\ 
—  =  d  .  (i  -  X  ) 
q  X  p' 


(78) 


for  K,  :  M  =  d 

1  y 


cn  ,2  ,  o  Oi 

^  d  (£y-  /p) 


^2  :  Ny  = 


_  yO) 

q  y  ‘'p/ 


where 


.1 


K,  :  0  =  d 
5  y 


.Ov2 


cn  j6  /  .0  o  \ 
—  =  d  (£  -  y  ) 

q  '  y 


-  x“)‘  +  (i“  -  y°)2  =  =  C 

Consequently  we  obtain; 


xy 


(for  the  meaning  of  C  and  C 

y  ^ 

see  formulas  (40)  and  (4l)) 


II 

o 

• 

C 

y 

1 

It 

•< 

O 

K 

=  . 

C 

2 

-N  =  C  .  C 

X 

xy 

y 

► 

>> 

=  • 

c 

-0  =  .  c 

x 

xy 

y 

y  xy  ,  ; 

(79) 
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The  least  squares  solution  gives,  besides  the  orientation  xonknowns,  the 
distortion  coefficients  K^,  and  K^.  The  distortion  A'  can  now  be  computed 
with  formula  (73)-  However,  usually  the  distortion  curve  is  presented  in  such 
a  way  that  A  =  0,  for  a  suitably  chosen  3.  Thus  analogous  to  formula  (71 ),  the 
condition  that 


K 


*Jl  .  +  Kf) 


(80) 


must  be  satisfied. 


From  (73)  the  right  side  of  formula  (80)  equals  -A* 

a 

and  the  final  distortion  curve  is  now 


*  5  5  7 

A  =  K^d  +  K^d-^  +  +  K  <!' 


(81) 


or 


r  *  2  4  6 

A  =  d  k  +  +  Kjd 


(82) 


The  "focal  length"  associated  with  this  distortion  curve  is  again,  in  analogous 
to  (71), 

(85) 


*  / 

c  -  c(l  -  K^) 


Accuracy  Considerations: 

The  mean  error  of  an  observation  of  unit  weight  denoted  by  m  is  computed 
according  to  (58) 

■  «  =  (8^) 
where  r  is  the  number  of  observation  equations 


u  is  the  number  of  unknown  geometrical  parameters  and 

Kis  the  number  of  unknown  distortion  parameters  carried  in  the  solution. 


As  described  in  Chapter  IV  E.  the  Inverse  of  the  normal  equation  matrix 
is  the  weight  matrix  of  the  unknown  parameters  of  the  solution  and  consequently 
the  mean  errors  of  the  distortion  parameters  can  be  computed  directly. 
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The  mean  errors  of  a  computed  distortion  A  for  an  arbitrarily  chosen  d 
remains  to  be  computed.  _ 

V  “  V  ^A  ( 

Substituting  (6l)  into  (82)  we  obtain 

A  =  d  +  K2(d^  -  +  K^(d^  - 


which  may  be  written  as 


denoting 


=  d  [f^K^ 


+  f^K^  +  f,K, 
2  2  5  5_ 


whereby 


we  may  write  (85)  as 


f^  =d(d2  - 
f^  =  d(d*^  - 
f^  -  d(d^  - 


r  1 

"A-”irV|frj 


where  Q|^  is  the  sub-matrix  of  the  Inverse  of  the  original  normal  equation 
system  associated  with  the  unknown  distortion  parameters  and  K^. 

VII.  CONSIDERATION  OF  REFRACTION- 

I  , 

Similar  to  distortion,  refraction  causes  the  tangent  to  the  actual  light 
ray,  at  the  center  of  projection,  to  deviate  from  a  line  expressing  the  con¬ 
dition  of  co- linearity  between  object  point,  center  of  projection  and  image 
point.  . 


In  order  to  compensate  for  refraction,  it  is  necessary  to  Incorporate  a 
corresponding  correction  in  the  computation  of  the  x°  and  y°  values  in  formulas 

/TON  +  ^  -ira1  'C'^  rrt-ivta,  *l  t*-«  +Vi.i  V>nn-4n  a-h  aaai.>uJ  aa 

^  WJ.  V  •  J.-  w  Juvy  DX1«^W(9  k^ciiOXW  ox  LlUA  a.pg  14111  L  J 


We  obtain  from  Fig. . (lO) ; 

ctn  z  - - - - —  (88) 

^  \/(X)^  +  {Yf 

A  denotes  the  correction  to  the  zenith  angle  due  to  refraction, 

A  =  F(X^  Z^,  8,;  z,  M)  (89) 

M  are  meteorological  parameters  and 

I  ■  ■  ' 

6  =/  [(X)^  +  (Y)^  +  (Z)^] 

Consequently,  we  obtain; 

> 

1. 

(Z)^  =  [(X)^  +  (Y)^]  ctn  (z  -  A)  (90) 

The  computation  of  x°  and  y°  with  formulas  (12)  is  now  carried  out,  using  for 
each  individual  ray,  the  corresponding  (x),  (Y),  and  (Z)^  coordinates.  As  the, 
(x)(y)(z)  values  converge  during  the  iteration  cycles  to  the  final  answer,  so 
will  the  corresponding  A  correction  converge  to  the  correct  refraction  value. 

Fig,  11  shows  the ^general  character  of  the  refraction  values  as  they  must 
be  expected  by  an  observer  bn  the  ground  (dotted  lines)  eind  in  the  air  (solid  . 
lines ) , 

Refraction  In  seconds  of  arc  Is  presented  In  Its  functional  relation  to 
the  elevation  or  depression  angle  of  the  line  of  sight  and  to  the  height  above 
the  reference  ellipsoid  of  target  point  or  observer,  respectively, 

“v-The  tables  No.  1  and  No.  2  show  the  same  Inforaatlon  In  somewhat  more 
detail.  The  columns  headed  by  A/1°  show  the  changes  of  refraction  for  1°  change 
of  a  specific  elevation  or  depression  angle,  respectively. 

It  appears  practical  to  assume  that  for  aerial  precision  measurements  the 
210nim  -  60°  l^ns  cone  under  20°  tilt  presents  the  most  stringent  requirements 
with  regard  to  refraction.  Therefore  55°  for  a  minimum  depression  angle  appears 
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to  "be  a  practical  limit.  Tatile  No.  2  shows  that  for  an  aerial  observer  both 
the  absolute  amount  of  refraction  and  its  differential  change  have  their  maximum 
values  for  a  flying  height  of  about  l8  km.  For  a  depression  angle  of  55°  the 
corresponding  values  are  27”  and  l",  respectively.  Based  on  practical  experi¬ 
ence  with  metric  aerial  photography,  one  should  be  satisfied  with  a  compu- 
tational  accuracy  of  +1"  for  the  refraction  correction.  Consequently,  the 
corresponding  depi^ssion  angles  for  the  Individual  rays  can  be  computed  by 
simple  formulas  depending  on  a  local  Cartesian  coordinate  system,  provided  that 
its  origin  is  within  +1°  =  +100  km  from  the  camera  position.  This  approach 

will  do  Justice  to  strips  of  200  km  length  or  blocks  with  sides  of  about  I50  km. 

At  the  same  time  it  is  desirable  to  minimize  the  computational  difficulties 
connected  with  the  introduction  of  certain  metric  conditions  associated  with 
partial  control  points.  A  corresponding  orientation  of  the  local  Cartesian 
system  has  been  described  in  Chapter  IV  -  (D),  page  55, 

Formulas  describing  refraction  corrections  have  been  derived  by  several 
authors.  In  Cs]  a  treatment  is  given  for  astronomical  refraction  including 
zenith  distances  ^  90°,  In  expressions  for  target  points  inside  the 
atmosphere  have  been  derived.  References  with  respect  to  corresponding  basic 
assumptions  and  formulas  may  be  found  in  the  recently  published  report,  [iC^  . 

A  summary  of  the  problem  of  refraction  in  photogramme try  is  found  in  Cli]  . 

M,  the  meteorological  constants  referred  to  in  formula  (89)  are  either 
obtained  in  accordance  with  some  model  of  the  atmosphere  or  are  determined 
directly  from  Independent  meteorological  measurements,  or  obtained  indirectly 
by  additional  elevation  angle  measurements  with  respect  to  known  points;  a 
procedinre  which  under  certain  conditions  can  he  made  a  part  of  the  actual 
photogrammetrlc  triangulation  measurements ,  If  aerial  analytical  triangulation 
demands  specific  meteorological  parameters  in  the  .area  of  the  flight,  it  is 
feasible  to  eject  periodically  from  the  airplane  probes  for  measuring  densities 
and/or  temperatures;  the  corresponding  results  being  recorded  by  radio  link  in 
the  airplane. 
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In  accordance  vith  the  conventional  approach  we  derive  from  Fig.  (13)  ■ 

and  Snell's  law: 


sin  (z)  n 
P  P 


n  =  index  of  refraction 


sin  ^ 

and  from  the  sine  law: 


(91) 


sin  (z)  r  ,  /  v 

P  _  P  where  r^  =  R  +  (Fig,  12) 

sin  t  r 


(92) 


with  (91)  a^d  (92)  it  follows  that; 


*  *  * 

n  r  sin  (z)  =  n  r  sin  (z)  =  K 

P  P  P  P  P  .  ^P 


(95) 


in  particular  for  point  A,  K  =  n  r  sin  (z)  where  r  =  R  +  H 
Further  from  Figure  (15),  omitting  subscripts. 


dC 


K  dr _ 

r  Qnr)^  - 


Applying  formula  (94)  to  points  (Q)  outside  of  the  effective  atmosphere,  where 

n  =  unity,  we  obtain  with  r  =  R  +  H 

1  <1 


a 

I 


(95) 


and,  correspond.ingly,  from  Fig.  (12)  for  all  points  for  which  r  S  r 

'  b 

C  -  +A00  - 

<1 


(96) 
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TECKKICALlU„.^^j 
BRANCH  NO.  1 

From  (96)  It  follovs  that,  for  target  points  situated  outside  the  effective 
atmosphere,  refraction  can  he  expressed  In  terms  of  astronomical  refraction 
denoted  hy^o^j.  According  to  [8]  we  have 

*  3  *  5  *  7  * 

,  tan  z  ’  +  tem'  z  +  tan^  z  +  tan'  z  +  . . .  (97) 


with  ctn  2z  =  V  ctn  (z)  and  for  air  v  =  Q.1^78\l^^ 
a  ''a  a  V^a 


The  coefficients  depend  upon  the  structure  of  the  atmosphere.  As 

outlined  In  for  precision  work  It  will  he  necessary  to  compute  with  the 

Cauchy  equation  as  function  of  the  effective  wave  length  and  obtain  an  Index;  of 
refraction  profile  from  Rawlnsonde  observations .  It  Is  then  possih'ie  with  . 
formula  (9^)  to  determine  for  varioiis  (z)  values  the  corresponding  gj- values  by. 
numerical  integration.  The  corresponding  /\  values ,  computed  from  the  re¬ 
lation  =  (z)  +  (Fig.  12)  are  now  used  to  compute  the  specific  a- 

coefficlents  by  fitting  the  expansion  as  given  by  (97)  to  the  computed 
values . 


In  order  to  apply  formiilas  (9^)  or  (96)  to  the  case  of  an  aerial  observer, 
it  Is  necessary  to  establish  the  relation  between  (2;)^  and  (z)^  which  Is 
obtained  from  (95). 


sin  (z) 


K 


n  r 
a  a 


n  r  sin  (z) 
p  p  '  P 

n  r 
a  a 


(98) 


The  true  zenith  distances  z  and  z 

a  p 


are  computed,  according  to  Fig.  12  from 


and 


tan  z 

a 


sin  C 

_ _ E 


cos  C 

P 


r 

a 

r 

P 


(99) 


'1 


tan  z 

P 


-sin  C 
_ 

r 

cos  C  -  — £ 


(100) 
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From  Fig.  i2  we  see  that  for  all  points  outside  the  effective  atmosphere^ 
the  refraction  as  encountered  hy  an  observer  on  the  ground  and  in  the  air  for 
a  specific  ray  add  up  to  astronomical  refraction. 

This  fact  is  illustrated  in  Fig.  11  where,  as  example,  for  an  elevation 
and  depression  angle  of  50°^  the  sums  of  the  ordinates  of  the  solid  and  dotted 
lines  have  been  plotted  for  altitudes  outside  of  the  effective  atmosphere.  The 
resulting  line,  designated  "astronomical  refraction  for  50°",  agrees  with  the 

corresponding  value  of  99"  in  Table  No.  1  for  c  =  50°  and  H  -  GO  .  i 

.... 

For  precision  triangulation  of  points  within  the  atmosphere,  especially  if 
large  distances  between  observer  and  target  are'  encountered  as  in  guided  missile 
applications,  the  numerical  integration  by  formula  (94)  based  on  an  index  of 

refraction  profile  will  be  unavoidable.  For  many  cases,  however,  especially  in 

< 

aerial  photogrammetry,  an  expression  for  the  refraction  can  be  used  which'  avoids 
the  cumbersome  nimoerical  integration. 

Assuming  a  constant  atmospheric  temperatiire  gradient,  refraction  ex¬ 
pressions  were  derived  in  [  9]  for  both  an  observer  on  the  ground  and  in  the  air. 
However,  these  formulas  are  restricted  to  a  situation  where  one  of  the  end 
points  of  the  light  curve  is  situated  at  the  height  of  the  ellipsoid  of  refer¬ 
ence.  -"Moreover,  the  expressions  obtained  inf  9l  unnecessarily  complicated 
for  numerical  evaluation.  In  the  following  a  derivation  is  given  'vrtilch,  follow¬ 
ing  the  general  approach  of  [9]  r  overcomes  the  above  mentioned  restriction  and 
results  in  expressions  more  suited  for  numerical  evaluation."' 

From  the  well  known  pressure  altitude  relationship  it  follows  that  for  a 
linear  decrease  of  temperature  with  altitude. 


T  \a 


£_  .  (Z_) 

p  'T  ' 

o  o 


where  a  -  »  ^  and  T  =  T  +  LH 
RL  o 


(101) 


and 


R  is  the  gas  constant  for  air,  T  absolute  temperature 
-dT 


L  = 


dH 


the  temperature  gradient. 


Similarly  well-lmown  is  the  density  -  altitude  -  pressure  relationship. 

f-  -  I-  r-  (-f)"' 


a-1 


Finally  we  introduce  the  well  known  relationship  between  the  indexes  of 

2 

refraction  and  density  (compare  e.g.  )  n  -  1  =  cp,  where  c  Is  a  constant. 

With  n*=  1  +  a,  where  a  Is  usually  <0.0003;  we  may  write; 


n-1  _ 

n  -1  ” 
o 


=  £- 


or  n  =  1 


/T  xa-l 


(105) 


( 


(104) 
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Prom  Pig.  l4  : 


X 


(105) 


=  (h  +Ah)  tan  (z) 


Ah  sin  (z)  Ah  sin  (z)  cos  (z) 

and  A  =  - -  w  - 5 5 

a  6  .  h 


(106) 


In  [»]  an  expression  for  x  is  published,'  which  may  be  written  as: 


tan  (z) 


H 


X  =  h  tan  (z)  +  „ 

'  'a  2 


COB  (z) 


/p  n  -  n 

- E-dH 


(107) 


E 


For  optical  frequencies  n^  Is  usually  <  1.0005  but>l  and  we  may  therefore 
write  for  (I07) 


X  =  h  tan  (z)  + 


tan  (z)^  ®p 


COB  (z) 


/ 


(n  -  n  )dH 

'  &  P' 


(108) 


By  comparing  (I05)  with  (108)  it  follows  from  (106) • 


tan  (z)^  JJp 

— —  J  K  -  V 


dH 


(109) 


Due  to  the  linear  relationship  between  temperature  and  height,  (lOl),  we  may 
write  with  the  notation  in  Figure  l4  and  formula  (l04').  : 


T 


tan  (z)  p 

— —r -  f  (n  -  n  )  dT 

t  j  '  a 


(no) 


From  (105); 


n  =  1  +  a  (^) 
a  o'T  ' 

o 


T 


n  =  1  +  a  (=^) 

P  o -T  ' 

o 


P\  .a-1 


_ 2 _  a-1  a-1 

(“a  -  "p)  -  ^-1  c-  ’■r")  =  ^r-  [1-  ]  nil) 
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Substituting  (ill)  into  (llO)  we  have: 


A  = 


tan  (z)  a  T 
a  o  a 


,a-l 


T 


a-1 


0 


1  -  (1  + 


dt 


A  = 


tan  (z)  a  T 
'  a  o  a 


,a-l 


„a-l 


.  a  t  ^a 

t - (1  +  — )  + 

a  '  T  ' 

_  a 


iA  „  =  — tan  ( z ) 
■^o 


T 


,a-l  a  p 


a(T  -T  ) 
^  a  p' 


or 


A  =  A  tan  (z) 
a  ^  'a 


where 


n  1 

_  a  p 


a(T  -T  ) 

a  p'  _J 


a 


A  = 


a  =  - 


o 

i_ 

RL 


T  =  T  +  LH 
a  o  a 


T  =  T  +  LH 

pop 


From  Figure  .12  or  l4-  we  have; 


T  =  A  +A 
a  p 


(112) 


(113) 


(114) 


(115) 


(116) 


where  t  is  angular  refraction. 

From  93  follows : 

n  r  sin  (z) 
a  p  '  ^  a 

n  ~  r  sin  (z) 
pa  '  ''p 

r 

with  (ll6),  setting  — 

a 


and  A  is  the  refraction  of  an  aerial  observer. 
P 


(117) 


1,  and  with  t  being  a  small  angle. 


we  may  write  (117)  s-s: 
n 

—  =  1  +  T  ctn  (z) 
n  '  '^a 

P 


(118) 
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(119) 


Neglecting  again  n^  in  tKe  denominator,  we  have: 

(n  -  n  )  =  T  ctn  (z) 

'  a  p'^  '  ''a 

or  with  (ill). 


O  ,  /  \  /  —Qr^l  _Sr"l\ 

■^o 


(120) 


and  with  the  notation  of  formula  (115) ^ 

t”  =  A  tan  (z)  (T®'“^  -  T^“^) 

'  ‘^a  '  a  P  ' 


(121) 


The  refraction  for  an  aerial  observer  follows,  from  (ll6): 


A"  =  A  tan  (z) 
p  'a 


T  -  T 

_ P_  _  rpa-l 

a(T  -T  )  P 

'  Q  -1^  ' 


(122) 


As  mentioned  on  page  64,  t  becomes  astronomical  refraction  (^gg)  for  all  points 

outside  the  effective  atmosphere  for  which  =  0.  Consequently  (I2l)  reduces 

for  T  =0  to; 

P 

A"  -  A  tan  (z)  (l25) 


Introducing  (125)  into  (ll5)  we  obtain; 


A^'  =A^  -  A  tan  (z)^ 


_a  a 

T  -  T 
a  p 

a(T^-Tp) 


(124) 


-  T®" 

a  P _ 

aT®"'^(T  -T  ) 
a  '  a 


(125) 


whereby  for  precision  workA^  can  be  computed  from  (97). 

Refraction  can  now  be  taken  into  consideration  by  computing  either  A  or 
( z  -  A )  as  required  In  formula  ( 90 ) , 
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'  ^ 

Type  of  application 

Location  of  observer  or  targe' 

Outside  effective  atmosphere 

b  point  respectively. 

Inside  effective  atmosphere 

Aerial  Photogrammetry 

Formulas  (lOO),  (96) 

(97)  and  (95) 

Formula  (122) 

Terrestrial 

Photogrammetry 

Formulas  (99) >  (9^), 

(97)  and  (95) 

Formula  (115)  or  (124) 
and  (97)  or  (99)  &  (94) 
in  order  of  Increasing 

precision 

Attention  must  be  given  to  the  fact  that  the  computation  of  z  for  a  specific 
ray  is  carried  out  with  the  coordinates  of  the  center  of  projection  and  of  a 
specific  control  point.  Thus  z,  and  not  (z),  is  being* obtained  (Fig.  12).  In 
extreme  cases  of  precision  terrestrial  photo gramme try,  it  may  therefore  become 
necessary  to  compute  a  series  of  refraction  corrections  according  to  the  follow¬ 
ing  steps: 

1.  A*  =  f  z 

2.  {z)'=  z  -  A ’ 

5.  A”  =  f(z)  *  etc. 

The  computations  are  continued  until  the  specific  A  -value  stabilizes . 

Provided  (z)  +  (z)  is  sufficiently  close  to  l8o'^  (flat  earth  geometry), 

3,  p 

formiolas  (115)  and  (122)  can  be  combined  in  an  expression  valid  for  both  a 
ground  based  or  airborne  camera.  Introducing  the  notation  of  this  report,  the 
center  of  projection  by  0  and  the  target  point  by  R,  we  obtain: 


A" 


A  tan  (z) 


T 


,a-l 


T^  -  T^ 
^0  -^R 


(126) 
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VIII.  THE  TREATMEMT  OF  STRIP  AMD  BLOCK  TRIANGULATION 

It  has  been  shown  that  the  matrix  of  normal  equations  (35),  illustrated  in 
Fig.  15,  expresses  the  general  photogrammetric  problem.  The  contents  of  that 
matrix  Include  the  solution  for  the  trlangulatlon  of  an  unlimited  block  and 
obviously,  as  a  special  case,  the  trlangulatlon  of  an  unlimited  strip. 

Such  matrices  resemble  each  other  In  the  geometric  arrangement  because 
successive  groups  of  ground  points  are  being  photographed  from  a  certain  number 
of  successive  camera  stations .  This  fact  Is  reflected  In  the  arrangement  of 
the  coefficients  In  the  bJ(AP  niatrlx  and  Its  transpose.  These 

matrices  are  only  partially  filled,  forming  an  escalator  pattern  as  shown  In 
Fig.  15. 

The  example  Is  a  strip  with  2/5  overlap,  where  It  Is  assumed  that  six 

points  are  located  In  each  trllap-area  on  the 'ground  (l,  II . N).  The 

corresponding  reduced  normal  equation  system  (37)  Itself  becomes  a  symmetrical 
escalator  matrix,  grouped  along  the  main  diagonal  as  shown  In  Fig.  I6. 

Fig/  17  shows  a  reduced  normal  equation  matrix  for  the  determination  of 
the  corresponding  vectors  of  a  block  of  7  x  7  photographs  flown  with  2/5 
longitudinal  and  2/5  lateral  overlap.  In  this  way,  each  portion  of  the  ground 
Is  being  recorded  on  nine  photographs.  In  the  exanqjle  It  was  assumed  that  one 
point  Is  located  In  each  one  of  the  nine  times  covered  ground  sections.  The 
matrix  as  shown  In  Fig.  17  Is  typical  In  Its  arrangement  for  any  block  trl¬ 
angulatlon. 

The-  size  of  the  cross  shaped  openings  which  are  filled  with  zeros  obvi¬ 
ously  increases  with  an  Increase  In  the  length  of  the  sides  of  the  block  under 
consideration. 

If  one  uses  the  aforementioned  schematic  of  the  general  solution,  present 
day  electronic  computers  allow,  without  undue  difficulties,  the  formation  of 
the  reduced  normal  equation  systems  as  presented  In  Figs.  16  and  17.  However, 
the  problem  remains  to  Invert  these  normal  equation  systems.  Generally 
speaking,  both  strip  and  block  trlangulatlon  will  result  In  normal  equation 
systems  with  too  many  unknowns  for  direct  Inversion. 
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IFIQURE  17 


Two  possibilities  for  the  solution  of  this  problem  shall  be  discussed. 

The  first  solution  is  rigorous.  Its  result  is  Influenced  only  by  the 
propagation  of  the  rounding  errors  in  the  computing  machine. 

The  solution  depends  on  the  repeated  application  of  formulas  (22)  -  (50) » 
th]ds  eliminating,  stepvlse,  certain  groups  of  unknowns.  The  system  of  peu:- 
tltlonlng  as  displayed  in  Fig.  l6,  Is  based  on  the  stepwise  elimination  of 
groups  containing  twelve  unknowns .  This  arrangement  was  chosen  in  order  to 
minimize  the  nuinber  of  zeros  in  the  individual  submatrices,  denoted  by  B^'s, 
thus  increasing  the  economy  of  the  necessary  computations.  Denoting  the 
s -times  reduced  system  of  normal  eq^uations  in  accordance  with  the  well  known 
notation  used  by  Gauss,  as  .s, 

we  obtain:  =  [  Ag  •  (S  - 1 )]  '  [^ 'CS- 1 )] 

[^-t “  ®S-I  ^Og] 

I  ,  • 

^0|  '  [-1^1  “  ^02] 

(127) 

where 

AjI  =  Aj-b'i  A‘'B|  and  Pjl  = /j  -  B',  A7' J} 

Aj  Z  =  Aj- B2(A2  ir' B2  Pj  Z  = 

As-(S-I)  -As-B;.,[As.,-(S-2)]''Bs.,  4(S-l)  =  4-B^.l[As.,(S-2f4.|<S-2) 
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The  other  solution  uses  a  technique  related  to  the  Gauss-Seidel  relaxation 
method.  A  submatrlx  chosen  for  geometrical  reasons,  jnoves  in  steps  across  the 
original  matrix.  Once  more  the  escalator  matrix  of  Fig.  l6  is  seen  in  Fig.  l8 
as  the  shaded  area  within  the  heavy  contour  lines.  In  a  strip  flown  with  2/3 
overlap,  five  consecutive  camera  stations  are  connected  by  the  resulting  overlap 
of  the  photographs.  Therefore,  having  six  unknowns  per  station,  a  30  x  50 
submatrlx  was  chosen,  which  always  contains  such  a  group  of  unknown  orientation 
parameters,  as  they  belong  to  five  consecutive  camera  stations.  This  partial 
system  is  now  displaced  along  the  strip  by.  one  station  each  time.  Thus,  for  a 
strip  with  n  photographs  one  obtains  (n-h)  such  submatrices  and  consequently 
by. inversion  of  these  submatrices  one  obtains  five  values  for  each  of  the 
orientation  parameters  with  the  exception  of  the  first  and  last  four  photographs, 
where  accordingly  fewer  values  are  obtained.  The  arithmetic  means  of  the  roots 
of  the  individual  parameters  are  now  computed  and  considered  as  the  result  of 
any  one  iteration  cycle.  The  approximation  results  thus  obtained  are  used  to 
continue  the  computation  according  to  the  Gauss-Seidel  relaxation  method,  by 
changing  the  original  absolute  column,  taking  into  account  the  coefficients  not 
incorporated  in  the  individual  50  x  50  submatrices,  together  with  the  values  of 
the  corresponding  orientation  parameters,  as  obtained  in  the  preceding  iteration 
cycle.  The  Iteration  is  continued  until  the  roots .have  converged  to  a  pre- 
established  accuracy  level.  For  the  economy  of  the  solution,  it  is  of  importance 
that  only  in  the  first  iteration  cycle  the  individual  submatrices  must  be  in¬ 
verted.  The  roots  in  the  following  cycles  are  then  determined  by  multiplying 
the  individual  inverses  by  the  changed  absolute  columns. 


IX.  SUMMARY  AM)  CONCLUSIONS 
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The  analytical  solution  for  the  general  problem  of  photogrammet'ry,  as 
presented  in  this  report,  is  not  restricted  by  geometrical  or  statistical 
considerations,  because  all  nine  geometrical  parameters  which  characterize  a 
central  perspective  can  be  introduced  for  jany  number  of  photographs.  (Compare 
[l6] ) .  Furthermore,  provisions  have  been  made  to  consider  all  types  of  measure¬ 
ments,  as  they  may  arise,  as  erroneous.  A  least  squares  treatment  results  in 
the  most  probable  values  of  the . unknowns  of  the  solution,  provided  that  the 
residual  errors  are  normally  distributed  and  the  varlo\is  bundles  of.  rays  are 
generated  according  to  the  principle  of  the  central  perspective. 


Physical  Influences,  which  deform,  the  central  perspective  bundles  and 
therefore  devaluate  the  condition  of  colllnearlty  (compare  Chapter  III),  are 
lens  distortion  and  refraction.  The  lens  distortion  can  he  determined  along 
with  the  geometrical  parameters  (See  Chapter  VI ).  Star  or  collimator  pho¬ 
tography,  in  connection  with  the  analytical  treatment  of  an  individual  photo¬ 
graph,  provides  a  practically  unhlased  means  of  determining  the  calibration 
of  such  a  cameraj  that  la,  its  interior  orientation  and  the  specific  lens 
distortion.  The  simultaneously  obtained  elements  of  exterior  orientation  allow 
the  calibration  of  phototheodolltes  or  the  establishment  of  the  mutual  relation 
between  several  cameras  as,  e.g.,  is  necessary  in  connection  with  the  use  of  a; 
sun  camera.  (See  Chapter  IV -C) 

Based  on  practical  experience,  it  can  be  generally  stated  that  the  analyti¬ 
cal  treatment  appears  to  be  an  excellent  means  for  analyzing  objectively  the 
con^ionents  of  the  photogrammetrlc  measuring  method. 

If  refraction  can  not  be  eliminated  by  a  suitable  arrangement  of  the 
meastirlng  set-up,  its  Influence  must  be  eliminated  for  each  individual  ray, 

ng  use  of  specific  meteorological  measurements  or  assuming  a  certain  model 
of  the  atmosphere.  Information  concerning  the  computation  of  the  corresponding 
corrections  their  consideration  in  the  Einalytical  reduction  method  is  given 
In  Chapter  VII. 

Ignoring  the  economical  side  of  the  problem,  in  precision  photogrammetry 
the  recmirements  for  flatness  and  dimensional  stability  of  the  emulsion  carrier 
can  be  satisfied  by  the  use  of  precision  ground  glass  plates.  The  unavoidable 
irregular  shrinkage  of  the  emulsion,  together  with  the  measuring  errors  on  the 
comparator  cause  the  plate  coordinate  meaaiirements  to  be  affected  by  residual 
errors,  which  have  a  distribution  similar  to  a  normal  distribution,  thus  justi¬ 
fying  the  effort  of  a  rigorous  least  squares  treatment.  Such  a  computing 
technique  will  also  prove  to  be  practical  in  cases  where  systematic  errors  are 
present.  Wot  only  does  a  least  squares  solution  in  any  case  produce  the  most 
likely  result,  but  addltloneJJLy  obtained  qualitative  and  quantitative  infor¬ 
mation  allows  one  to  recognize  and  Isolate  systematic  errors.  In  addition, 
the  least  squares  treatment,  as  described  In  this  report,  is  simple,  due  to 
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the  simplicity  of  the  mathematical  model  on  which  the  general  photograrametrlc 
solution  is  based.  Thus,  the  least  squares  solution,  can  he  considered,  an 
economical  technique  for  intersecting  corresponding  rays  of  various  bundles. 

A  comparison  of  this  report  with  [  2]  ,  shows  that  the  present  approach 
distinguishes  itself  from  the  former  one  only  by  the  way  in  which  the  ellml- 
nation  of  the  unknown  coordinates  of  the  model  X  is  accomplished.  In  the 
present  solution,  this  elimination  of  unknowns  Is  not  performed  algebraically 
before  the  observation  equations  are  formed  but  during  the  process  of  forming 
the  reduced  normal  equations.  Consequently,  the  system  of  reduced  normal 
equations  (formula  (57))  agrees,  to  all  but  second  order  terms,  which  have  been 
neglected  during  the  process  of  linearization  by  the  Taylor  series,  with  the 
corresponding  normal  equation  system  of  the  formerly  published  solution.  The 
advantage  of  the  present  solution  can  he  seen  in  the  simple  and  systematic  flow 
of  the  computations.  It  is  possihS^  to  treat  all  cases  of  practical  analytical 
photogrammetry  with  hut  one  basic  computing  scheme,  thus  simplifying  consider¬ 
ably  the  "bookkeeping  effort"  in  the  electronic  computer, 

A  critical  study  of  the  indivual  steps  of  the  presented  solution  leads 
to  a  conclusion  which,  although  somewhat  discouraging  for  the  author,  may  en¬ 
courage  the  application  of  analytical  photogrammetry.  It  hecpmes  obvious  that 
the  analytical  treatment  of  photogrammetrlc  problems  does  not  call  for'  any  new 
manipulations  in  photogrammetrlc  theory  or  statistical  treatment  of  errors. 

The  expressions,  representing  the  basis  of  the  whole  solution,  derived  in 
equations  (ll)  and  (12),  as  they  exist  between  the  coordinates  of  the  model  anf} 
the  corresponding  coordinates  of  the  Images,  are  the  well-known  formulas  derived  , 
by  V.  Gruber  in  [Ik]  and  the  corresponding  inverse  functions.  The  partial  differ 
entlal  quotients  necessary  to  form  the  observation  equations  (15)^  given  wlth'< 
formulas  (40)  and  (kl)  are  the  same  expressions  that  are  found  in  [11  or  [21  ; 
the  unmodified  use  of  these  expressions  appears  justified  due  to  the  simplicity 
of  their  construction.  The  system  of  normal  equations  (2l)  resulting  from  the 
system  of  observational  equations  (l8)  is  identical,  as  mentioned  before,  with 
the  solution  given  by  Helment  In  [41  .  The  elimination  of  the  vectors,  K  In 
the  system  of  normal  equations  (2l),  In  the  system  (55),  and  the  reduction' 
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of  the  normal  equation  system  for  an  unlimited  strip  as  shown  in  Pig.  (l6), 
(formula  (127))^  are  performed  with  a  certain  sequence  of  matrix  operations, 
formulas  (22)  to  (30)..  This  technique  however,  is  as  shown  hy  Gotthardt  in 
[15] none  other  than  the  classic  Gausslon  elimination  for  two  unknowns, 
applied  to  matrix  calculus.  The  alternate  solution  for  solving  the  system 
of  normal  equations  of  a  strip  or  block  triangulation  using  an  iterative 
approach,  is  based  on  the  Gauss-Seldel  Relaxation  method [ij] ,  in  connection 
with  a  procedure  which  has  proven  useful  in  similar  applications,  known  as 
"Smoothing  with  Moving  Arcs." 

Consequently,  the  knowledge  of  the  classical  geometrical  considerations 
dealing  with  a  central  projection,  and  a  certain  familiarity  with  the  method 
of  least  squares  adjustment,  as  applied  in  geodesy,  suffice  to  solve  the  ana¬ 
lytical  problems  arising  from  the  application  of  photograrametric  measuring 
systems . 

As  previously  mentioned,  the  system  of  normal  e^iuations  (35)  schematically 
shown  in  Fig.  4,  is  typical,  for  the  most  general,  problem  of  analytical,  photo- 
grammetry.  The  system  is  readily  accumulated  because  the  corresponding  obser¬ 
vational  equations  are  based  on  the  simplest  geometrical  model  conceivable, 
expressing  for  any  type  of  control  point  the  condition  of  colllnearity  between 
control  point,  center  of  projection  and  image  point.  The  problem  of  arriving 
at  the  most  economical  and  feasible  method  of  reduction  in  analytical,  photo- 
grammetry  is,  therefore,  concerned  with  the  process  of  determining  the  roots 
of  this  normal  equation  system.  The  solution  presented  in  this  report  is 
based  on  the  formation  of  a  system  of  reduced  normal  equations  (formulas  (3?) 
or  (38))  by  a  rigorous  mathematical  method.  The  attractiveness  of  the  solution 
arises  from  the  presence  of  a  series  of  separated  square  matrices  along  the 
main  diagonal.  This  featxjre,  however,  is  lost  in  the  resulting  system  of 
reduced  normal  equations,  where,  quite  obviously,  all  the  remaining  unknowns 
are  more  or  less  correlated,  depending  on  the  geometrical  arrangement  of  the 
cameras.  Any  attempt  to  further  simplify  the  process  of  determining  the  roots 
of  the  normal  equation  system  (35)  as  it  is  readily  seen  from  Fig.  4,  miist  try 
to  preserve  all  of  the  fully  separated  square  matrices  along  the  main  diagonal 
during  a  complete  computational  cycle.  Such  a  result  la  obtained  if  a  co^ipu- 
tational  cycle,  using  the  relaxation  techniques,  is  established  around  the  point 
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which  separates  the  unknown  parameter  corrections  associated  with  the  model 
from  those  associated  with  the  camera  orientations  Aq  .  Assuming, 
as  a  first  step,  Aq  a  null  vector,  one  obtains  a  Aj^  vector  by  a  series 
of  inversions  of  5x3  matrices,  each  of  which  contains  the  coordinate  cor¬ 
rections  of  a  single  control  point.  Multiplying  this  Av  vector  with  the 
T  T  "I  ^ 

Bubmatrlx^  ^0^^^  ^  ,  and  adding  this  result  to  the  absolute  column,. 

Bi(AP  'A')'i  ,  a  Aq  vector  can  be  computed  by  inversions  of  a  series  of 
maximum  (9^9)  matrices,  each  of  which  contains  the  corrections  to  the  orien¬ 
tation  elements  of  a  single  camera.  This  new  Aq  vector  is  multiplied  with 
the  submatrix,  BJ(AP  and  added  to  the  absolute  column,  Bj(AP''AVi 

and  a  new  Aj^  vector  is  computed  using  the  aforementioned  computational  steps, 
which  in  turn  lead  to  the  computation  of  a  new  Aq  vector.  This  method  will 
converge  although  very  reluctantly.  A  geometrical  analogue  of  such  a  solution, 
although  not  entirely  descriptive,  leads  to  the  following  approach.  Starting 
with  certain  approximations  for  the  orientation  parameters,  sets  of  coordinates 
of  all  points  of  the  model  are  computed  with  formulas  (56)  and  (37)^  which  may 
be  subjected  to  an  after-treatment  according  to  formula  (39)  with  Aq  as  null 
vector.'  In  any  case,  the  maximum  size  of  the  matrix  which  must  be  inverted  is 
(3x3).  With  the  thus  obtained  spatial  coordinates  of  all  points  of  the  model, 
a  series  of  resections  in  apace  is  computed,  thus  obtaining  the  orientations 
of  all  cameras.  In  these  computations,  matrices  of  (9x9)  maximum  size  must  be 
inverted.  With  the  thus  computed  orientation  parameters,  a  new  set  of  coordi¬ 
nates  of  the  model  are  computed,  on  which  a  new  set  of  camera  orientations  can 
be  based.  By  repeating  these  two  phases,  alternately,  the  final  orientations 
and  correspondingly,  the  final  coordinates  of  all  raeasxired  points  of  the  model 
can  be  computed.  Again, ’ the  convergence  is  extremely  slow.  An  Increase  in 
the  slope  of  convergence  woiiLd  remedy  this  situation.  The  associated  numerical 
effort  may  be  considerable.  Such  methods  are  described  in  [.17].  On  the  other 
hand  the  development  of  electronic  computers  progresses  at  an  Impressive  rate, 
with  respect  to  both  storage  facilities  and  computing  speeds.  In  the  near 
future,  it  should  be  possible  to  handle,  economically,  numerical  solutions 
requiring  a  very  large  number  of  iteration  cycles ,  It  is  believed  that  this  ' 
situation  will  make  possible  a  solution  based  on  relaxation  techniques  for  the 
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general  photogrammetrlc  problem.  This  remark  is  made  especially  with  respect 
to  the  analytical  treatment  of  an  extended  block,  where  the  approach  as  outlined 
in  this  report  still  leads  to  a  rather  bulky  matrix  of  normal  equations,  as 
shown  in  Fig.  17. 

Concluding,  it  appears  that  the  problem  of  analytical  photogrammetry  today 
requires  concentration  not  so  much  on  the  problem  of  the  numerical  treatment  of 
the  measured  plate  coordinates,  but  more  on  the  technical  difficulties  associ¬ 
ated  with  problems  of  Identifying  and  measuring  precisely  the  contents  of 
photographs . 

HELLMUT  H,  SCHMID 
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